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ABSTRACT 


Markov  Chain  models  have  been  used  to  forecast  stocks  in 
a  wide  range  of  manpower  systems.  Studies  have  been  done  in 
many  areas  such  as  education  planning*  hospital  planning, 
manufacturing,  private  research  and  development,  a  women’s 
military  unit,  the  civilian  work  force  supporting  the  U.S. 
Navy  and  a  state  police  organization.  This  study  looks  at 
such  systems  under  conditions  of  change  and  develops  the 
eguations  that  describe  the  steady  state  distribution  of 
personnel.  The  conditions  of  change  include  systems  where 
recruitment  is  constant,  increasing  (decreasing)  additively, 
or  increasing  (decreasing)  multiplicatively  and  systems 
where  the  changes  in  total  system  size  are  additively  or 
multiplicatively  increasing  (decreasing) .  Numerical  exam¬ 
ples  utilizing  these  models  are  provided,  along  with  a 
computer  program  of  the  formulas  written  in  the  language 
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I.  BACRG80PMD 


A.  INTRODUCTION 

According  to  Bartholomew  and  Forbes  [Ref.  1:  p.  1% 

Manpower  planning  is  often  defined  as  the  attempt  to 
match  the  supply  of  people  with  the  jobs  available  for 
them. 

This  is  a  concept  that  is  an  integral  part  of  military  force 
planning,  today.  All  branches  of  the  Department  of  Defense 
are  involved  with  studies  concerning  manpower  planning. 
Since  large  numbers  of  personnel  are  involved,  the  use  of 
manpower  models  to  determine  the  results  of  changes  to  the 
existing  system  before  they  are  actually  implemented  can 
prevent  costly  future  problems.  This  is  especially  true  in 
terms  of  having  the  right  number  of  personnel  in  the  rignt 
places  with  the  right  grade  levels  ana  skills,  at  the  right 
time.  To  quote  again  from  Bartholomew  and  Forbes  [Ref.  1: 
p.  6 1# 

A  manpower  model  is  a  mathematical  description  of  how 
change  takes  place  in  the  system. 

Much  work  has  been  done  in  the  area  of  manpower/ 
organizational  modeling.  The  text  by  Bartholomew  and  Forbes 
(referenced  above  and  hereafter  referred  to  as  BaF)  compiles 
information  cn  the  subject  from  numerous  disciplines  and 
standardizes  notation  to  describe  some  of  the  basic  models 
and  statistical  techniques  used  in  manpower  planning.  A 
central  theme  of  this  work  is  the  idea  of  depicting  an 
organization  as  a  system  of  stocks  and  flows  [Ref.  1:  pp. 
3,4].  The  personnel  within  the  organization  under  study  are 


divided  into  classes,  or  categories,  by  specific  character¬ 
istics  (for  example:  age,  time  in  service,  grade,  etc). 

The  numbers  of  personnel  within  categories  at  a  specific 
time  are  considered  stocks.  Flows  are  the  movements  of 
personnel  from  one  category  to  another  during  a  unit 
interval  of  time,  such  as  a  year. 

Various  mathematical  models  that  predict  future  stocks 
from  one  year  (or  other  time  period)  to  the  next  are 
described  in  BaF.  Of  particular  interest  here  is  the  tran¬ 
sition  model  that  is  based  on  the  theory  of  Harkov  Chains. 
In  addition  to  describing  the  model  [Bef.  Is  pp.  85-98],  BaF 
provides  a  BASIC  computer  program  [Bef.  1:  pp.  248-260  ],  to 
calculate  future  stock  vectors.  That  program  was  converted 
into  an  APL  program  at  the  Naval  Postgraduate  School  (NPS) 
and  is  available  on  the  NPS  mainframe  computer,  an  I3H 
system  370.  The  APL  program,  like  its  predecessor,  computes 
stocks  from  one  year  to  the  next  based  on  input  such  as, 
initial  stocks,  recruitment  and  wastage  (or  attrition)  rates 
and  the  transition  probabilities  for  '  the  stocks.  This  can 
be  done  repetitively  and  theoretically  could  compute  future 
stocks  forever. 

There  are  two  versions  of  the  basic  llarkov  chain  model 

that  are  well  developed  in  BaF.  One  is  based  on  a  fixed 

flow  of  recruits  into  the  system  and  the  other  on  a  fixed 
total  system  size.  Also  of  interest  are  other  versions  (or 
options)  that  deal  with  a  system  undergoing  growth  where 
either  recruitment  or  total  system  size  grows  in  an  additive 

or  multiplicative  fashion.  The  APL  program  allows  the 

computation  of  future  stocks  for  the  cases  where: 

(1)  recruitment  is  fixed,  grows  additively  or  multi- 
plicatively; 

(2)  the  total  system  size  grows  additively  or  multi- 
plicatively. 


Steady  state  is  the  condition  which  a  system  reaches, 
after  sufficient  time  has  passed,  so  that  the  stocks  no 
longer  change.  There  have  been  steady  state  models  devel¬ 
oped  for  the  cases  where  there  is  a  fixed  flow  of  recruits 
into  the  system  and  where  the  total  system  size  is  fixed. 
In  both  these  cases  the  number  of  personnel  in  each  category 
remains  fixed  after  steady  state  is  reached  by  the  system. 
The  vector  whose  components  are  these  numbers  is  called  the 
steady  state  stock  vector  (S.S.S.V. ). 

Under  conditions  of  growth,  the  steady  state  of  a  system 
must  be  redefined.  In  those  cases,  steady  state  does  not 
exist  in  the  above  sense,  because  the  system  is  undergoing 
constant  expansion  and  so  stocks  will  never  reach  a  fixed 
size.  However,  it  turns  out  that  the  stock,  sizes  do  reach 
constant  proportions  and  so  it  is  possible  to  talk  about  a 
steady  state  distribution  vector  (S.S.D.V.).  Components  of 
this  vector  represent  the  proportions  of  personnel  that 
remain  in  each  cateqory  forever  once  steady  state  has  been 
reached.  In  effect,  in  taose  cases  where  the  distribution 
of  personnel  among  the  categories  no  longer  changes  with 


time  the  system  has 
distrib  ution.  " 


reached 


steady  state  "in 


Althougn,  at  the  present  time,  the  APL  program  at  the 
Naval  postgraduate  School  predicts  future  stocks  for  a 
system  that  grows  under  conditions  mentioned  above,  it  is 
able  to  supply  a  steady  state  stock  vector  (S.S.S.V.)  only 
for  the  case  where  there  is  a  fixed  flow  of  recruits  into 
the  system.  The  only  way  to  get  the  steady  state  distribu¬ 
tion  vector  under  the  various  conditions  of  growth  is  to 
continue  to  run  the  program,  i.e.,  to  forecast  stocks  for  as 
many  years  into  the  future  as  necessary  until  steady  state 
is  reached.  This  is  not  a  satisfactory  method  for  several 
reasons  : 
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<1)  it  takes  too  long  to  run  the  model  when  the  system 
reaches  steady  state  only  after  many  years  have 
passed, 

(2)  it  is  a  trial  and  error  process  of  going  back  and 
forth,  in  time,  to  try  and  find  out  at  what  point 
in  time  steady  state  has  been  reached, 

(3)  occasionally,  even  if  the  program  is  run  out  as 
far  as  the  present  formatting  of  the  program 
allows,  there  can  still  be  instances  where  steady 
state  has  not  yet  been  reached. 

This  thesis  uses  the  *  basic  prediction  equation* ,  as 
described  in  BaF  [Hef.  1;  pp.  86-88  ],  to  correct  the  afore¬ 
mentioned  deficiencies  by  developing  models  to  compute 
analytically  the  steady  state  distributions  for  systems 
under  varying  conditions  of  growth.  Also  developed  was  an 
APL  program,  to  ne  added  as  a  subroutine  to  the  main 
program,  that  computes  the  S.S.D.V.  for  these  options. 

The  next  section  in  this  chapter  sets  out  the  necessary 
notation  and  equations  to  enable  the  reader  to  follow  the 
development  of  these  steady  state  models.  For  more  amplifi¬ 
cation  on  the  casic  Markov  model  see  BaF  [Hef.  1:  pp  3-8, 
86-132]. 

B.  NOTATION  AND  BASIC  EQUATIONS 

Stocks,  representing  personnel  at  an  organization,  are 
divided  into  k  categories.  The  stock  vector 


n(t)  =  (  n,  (t)  ,r. 2  (t) 


-  #  nK (t ) ) 


is  a  row  vector  where  each  element  n  •  (t)  ,  represents  the 
number  of  personnel  in  category  i  at  time  t.  Time  is  meas¬ 
ured  discretely  with  a  time  interval  of  unit  length  denoted 
(t-1,t)  .  The  most  frequently  used  period  of  time  is  the 
year  although  it  is  sometimes  convenient,  or  necessary,  to 
use  quarters,  months,  etc. 


Time  is  meas- 


Flows  describe  the  movements  of  personnel.  There  are 
two  basic  types  of  flows  in  a  system.  The  first  type  of 
flow  is  internal  to  the  system  and  can  be  thought  of  as 
transfers  from  one  category  to  another  within  the  system. 
Demotions  and  promotions  would  be  examples  of  this  type  of 
flow.  These  are  denoted  by  the  k  by  k  matrix,  N(t-1),  where 
each  element,  n(*j  (t-1),  represents  the  number  of  personnel 
moving  from  category  i  to  category  j  during  the  time 
interval  (t-1,t). 

The  other  type  of  flow  involves  the  two-way  transfer 
between  the  system  (organization)  and  the  outside  world. 
First,  there  are  wastage  (or  attrition  )  flows  of  the  number 
of  personnel  leaving  the  system  from  category  i,  during  the 
time  interval  (t-1,t),  denoted  by  n^  ^(t-1)  ,  (i  =  1,2,...  ,k 
and  t  =  1,2,...).  Then  there  are  recruitment  (or  accession) 
flows  of  personnel  into  category  j  from  the  outside  world, 
during  the  period  (t-1,t)  ,  denoted  by  n0j  (t)  ,  (j  =  1,2,...,* 
and  t  =  1,2,...).  Following  BaF  [Ref.  1:  pp.  3,4  ],  the 
notation  for  recruitment  flows  is  slightly  different  from 
that  for  internal  flows  and  wastage  flows.  This  is  to 

remind  the  reader  that  in  most  organizations  attrition  and 
internal  flows  (transfers)  are  accounted  for  first  and  then 
recruitment  is  determined  at  the  end  of  the  time  period. 

The  Markov  Chain  model  uses  flow  rates  (or  probabili¬ 
ties)  rather  than  actual  numbers.  P,  a  k  by  k  matrix  of 
probabilities  has  elements 


p‘j  * 


nj,j  (t-1) 
n l  (t-1) 


l ,  ^  1,2,... ,k  . 


There  are  certain  basic  assumptions  that  a  Markov  Chain 
model  needs  tc  meet  to  be  valid: 

(1)  the  transition  probabilities,  p^.  ,  do  not  change 

over  time;  ^ 


(2)  the  population  in  each  category  is  homogeneous  so 
that  p^«  represents  the  probability  of  each  indi¬ 
vidual  in  category  i  moving,  independently  of  any 
other  individual,  to  category  j. 

While  these  assumptions  may  not  always  be  met,  there 
have  been  enough  studies  done  using  the  model  to  show  they 
are  sufficiently  realistic  to  justify  the  use  of  the  model. 
Even  in  cases  where  these  assumptions  don*t  hold  very  well, 
the  model  was  fcund  to  be  quite  useful  for  predicting  future 
stocks  [Refs.  2,3]. 

Next,  w  =  (Wj  ,w2  ,...,wK)  is  defined  as  the  wastage 
(attrition)  rate  vector  where  each  element,  w;  ,  is  defined 


wi.  * 


n;  (t-1) 


Note  that# 

K 


2p‘j  *  1  • 


i  *•  1,2,..  .,k  . 


Finally,  r  =  (r,  ,r2  , -.  .  ,  rK )  is  defined  as  the  recruit¬ 
ment  proporticn  vector  where  each  element 


rJ  = 


nOJ  <*) 


is  the  proportion  of  total  recruitment  that  will  be  allo¬ 
cated  to  category  j,  with 

K  ’ 


=  2  S 


being  the  total  recruitment  during  the  time  interval 

C t—  1 , t )  .  This  implies  ^  rj  =  1. 

j=‘ 

BaF  [Ref.  Is  pp.  87,88],  snow  that  future  stocks  may  be 
calculated  using  the  eguation 


r 


.  (eqn  1.1) 


K 

nj  <*>  =  (t-1)p  * 

t*  i  J 

for  j=  1,2, ...,k  and  t  =  1,2,...  .  This  equation  states, 

that  the  stocks  in  category  j  at  time  t  are  equal  to  the 
number  of  personnel  that  moved  into  category  j  from  anywhere 
within  the  system  during  time  period  (t-1,t),  plus  the 
number  of  personnel  recruited  into  category  j  during  the 
same  period.  In  matrix  form  the  equation  is 

n(t)  =  n(t-1)P  +  R(t)r  t  =  1,2,...  .  (eqn  1.2) 

BaF  (Hef.  1,s  p.  88]  refers  to  this  as  "the  basic  prediction 
equation." 

Another  equation  dealing  with  changes  in  total  system 
size  is  also  developed  in  BaF  [Bef.  1:  pp.  94-96  ].  To  see 
this,  note  that 

K 

N(t)  =  2nc{t) 
i-  I 

is  the  total  system  size  at  time  t  and 

H(t)  =  N  ( t )  -  N  (t—  1 )  (eqn  1.3) 

represents  the  change  in  the  total  system  size  during  the 
time  interval  (t-1,t).  This  implies  that  total  recruitment 
during  (t-1,t)  can  be  expressed  as: 

K 

R  (t)  =  ((Wl»  ♦  «  (t)  *  0  -  (eqn  1.4) 

i*( 

The  first  term  in  this  expression  is  recruitment  that  is 
done  to  replace  those  who  left  (attrition  from  the  various 
categories)  while  the  second  term  consists  of  recruitment 
that  is  due  to  the  change  in  total  system  size. 

If  M (t)  =  0,  then  only  those  personnel  who  leave  the 

system  are  replaced,  which  implies  that  the  total  system 
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size  remains  fixed.  If  M  (t)  >  0,  new  personnel  are  supplied 
to  cover  an  increase  in  total  system  size,  i.e. ,  to  fill  new 
jobs.  If  M(t)  <  0,  then  some  jots  must  be  eliminated  as 
soon  as  they  are  vacated.  Of  course,  E(t)  >  0  always. 

Referring  back  to  Equation  1.1,  replace  R(t)  using 
Eguation  1.  4  and  let 

q*«  =  p  •.  ♦  w.  r  . 

‘J  ‘J  L  J 

Then,  as  BaF  [Ref.  1:  pp.  94,95],  shows 


nj  (t)  =  J  n<  (t_1)q//  +  M  <*>*; 

i *■/  * 

for  j=  1,2,...,k  and  t  =  1,2,... 


(eqn  1.5) 


To  put  this  into  matrix 


notation  let 


Q  =  (g  ■  • }  ,  i, j  =  1,2,. ..,k  , 
LJ 


Q  .  =  P  ♦  w  'r  . 


Then  Eguation  1.5  may  be  expressed  as 


n.(t>  =  _n(t-1)Q  ♦  M<t)£  ,  t  =  1,2, 


(eqn  1.6) 


Equations  1.2  and  1.6  form  the  basis  from  which  five 
steady  state  submodels  (or  options)  will  be  derived.  The 
next  section  introduces  these  options  and  the  direction  of 
development . 

C.  OBJ  ECTI YES 


The  objective  of  this  thesis  is  tc  develop  steady  state 
equations  modeling  five  different,  conditions  of  growth. 
Each  of  these  conditions  will  be  addressed  separately  as  a 
"growth  option."  The  first  three  options  deal  with 
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different  types  of  recruitment.  They  use  Equation  1.2  where 
R  (t)  is  the  total  recruitment  during  the  time  period  (t-1,t) 
and  in  these  three  options,  R  (t)  is  the  only  term  changing. 
The  last  two  options  deal  with  changes  in  total  system  size. 
They  utilize  Equation  1.6  where  M(t)  represents  the  change 
in  total  system  size  during  the  period  (t-1,t)  and  in  these 
two  options,  M (t)  is  the  only  term  that  changes.  A  complete 
list'  of  these  five  options  is: 

(1)  Fixed  recruitment:  fi  (t)  =  P,  t>1  ,  R>0  fixed, 

i.e.,  recruitment  is  always  constant; 

(2)  Additive  increase  (decrease)  of  recruitment: 

R(t)  =  R  ♦  (t-1)  M  ,  t>1 

where  a  is  the  constant  amount  of  yearly  increase 
(decrease)  ; 

(3)  aultiplicative  increase  (decrease)  of  recruitment: 

R ( t )  =  R0*-1,  t>1, 

where  100(0-1)  is  the  percent  of  yearly  increase 
(decrease)  ; 

(4)  Additive  increase  (decrease)  in  system  size: 

M(t)  =  a  ,  t>l 

where  a  is  the  constant  amount  of  yearly  increase 
(decrease)  in  total  system  size; 

(5)  aultiplicative  increase  (decrease)  in  system  size: 

M  (t)  =  (0-1)0*-»N(O)  ,  t>1 

where  100(0-1)  is  the  percent  of  yearly  increase 
(decrease)  in  total  system  size. 

If  a  =  0  in  option  2  or  0  =  1  in  option  3,  then  those 
options  revert  to  option  1,  of  fixed  recruitment.  If  a  =  0 
in  option  4  or  0  =  1  in  option  5,  then  fl(t)  =  0  and  the 
result  is  a  system  with  a  fixed  total  system  size.  That 
case  will  be  dealt  with  as  subcases  of  options  4  and  5.  The 
fixed  total  system  size  case  and  option  one,  where  R  is 
fixed,  are  the  only  cases  where  it  is  possible  to  achieve  a 


change  with  time.  In  all  other  cases  steady  state  will  be 
considered  to  have  been  reached  when  a  steady  state  distri¬ 
bution  of  personnel  is  achieved. 

Chapter  II  addresses,  in  depth,  the  analytical  develop¬ 
ment  of  the  steady  state  equation  for  each  of  the  options. 
It  starts  with  an  account  of  the  work  done  on  option  1  in 
BaF.  Then  it  proceeds  with  the  derivation  of  the  steady 
state  equations  for  options  2  through  5  from  either  Equation 
1.2  or  Equation  1.6  .  It  also  presents  examples,  one  a 
four-grade  hierarchical  system,  the  other  a  three  grade 
non-hierarchical  system.  In  addition  to  formulating  analyt¬ 
ical  models  for  options  2  through  5,  an  APL  program  has  also 
been  developed  that  will  compute  the  steady  state  distribu¬ 
tion  tor  those  options.  This  program  is  contained  in 
Appendix  F.  Appendices  A  through  E  are  computer  printouts 
of  terminal  sessions  of  the  two  examples  for  each  of  options 
1  through  5. 


II.  0PII08  DEVELOPMENT 


The  options  mentioned  in  Chapter  I  will  be  discussed 
here  under  the  usual  restrictions  of  manpower  systems. 
These  restrictions  imply  that  wj  >  0  for  all  i=1,2,...,k, 

since  in  manpower  systems  attrition  must  be  allowed  from  any 
category.  This  means  that  the  p  matrix,  as  defined  in 
Chapter  I,  is  composed  of  non-negative  elements  and  has  row 
sums  all  strictly  less  than  one.  The  Q  matrix,  also  defined 
in  Chapter  I,  is  a  stochastic  matrix  (i.e. ,  it  is  composed 
of  non-negative  elements  and  has  row  sums  all  egual  to  one). 
These  facts  along  with  the  following  two  mathematical  theo¬ 
rems  are  needed  in  the  derivation  of  the  steady  state 
results  for  the  five  options. 

Theorem  J 

If  P  is  a  matrix  composed  of  non-negative  elements  and 
row  sums  strictly  less  than  one  and  e  is  a  scalar  then 
the  matrix  0I-P  has  a  unigue  inverse  for  all  values  of 
0  >  1  [Bef.  4:  p.  42  ]. 

Theorem  2 

If  Q  is  a  stochastic  matrix  and  e  is  a  scalar  then  the 
matrix  ©I-Q  has  a  unigue  inverse  for  all  values  of 
0  >  1  [  Bef.  4:  pp.  6  0,6  1  ]. 

A.  RECBUITBENT  OPTIONS 

The  following  three  options  are  all  variations  of 
Eguation  1.2,  where  E  (t)  ,  representing  total  recruitment,  is 
the  term  that  takes  on  different  forms.  Conseguently,  the 
steady  state  equations  for  these  options  are  also  obtained 


in  somewhat  different  forms.  These  stead;  state  equations 
are  derived  in  the  next  three  subsections. 


1-  Potion  J:  Fixed  Recruitment 

The  steady  state  equation  for  this  first  option  is 
developed  in  BaF  [Ref.  1:  p.90].  The  analysis  presented 

here  is  a  brief  summary  of  that  in  BaF,  using  the  same 
assumptions.  Recalling  the  basic  prediction  equation, 
(Equation  1.2),  and  assuming  that 

limit  R  (t)  =  R  is  fixed 
t  —  oo 

it  follows  that  the  steady  state  stock  vector, 

n  =  limit  n  (t) 
t  —  oo  ~ 

also  exists  [Ref.  U:  pp.  40-43,48-50  ].  Therefore,  letting 
t—  oo  in  Equation  1.2  results  in 

n,  =  nP  ♦  Rr  , 


or 


n.(I  -  P)  =  Rr  . 

Finally,  knowing  the  matrix  I-p  has  a  unique  inverse  (see 
Theorem  1),  the  result 

n.  =  Rr(I  -  P)_l  ,  (eqn  2.1) 

is  obtained  for  tne  steady  state  stock  vector  (S.S.S.V. ). 
After  steady  state  has  been  reached,  it  follows  that 

n(t)  =  n_  , 

for  "large'*  t.  The  following  two  examples  will  be  used  to 
illustrate  the  use  of  this  steady  state  equation. 


This  first  example  uses  data  from  Example  4.6  in  3aF 
[fief.  1:  p.  97],  which  is  taken  from  the  women  officer's 
system  of  on^  of  the  British  services.  It  is  a  hierarch¬ 
ical,  four  category  system. 

The  inital  stock  vector  is: 

n(0)  =  (129,  74,  28,  11) 


The  transition  matrix  is: 


.  728 


.867  .033 


The  total  recruitment  is: 


H(t)  =  35  ,  fixed  for  all  t>1  . 
The  recruitment  proportion  vector  is: 


r  *  (1,  0,  0,  0)  . 


Then 


I  -  P  = 


-.046 


-.033 


and  the  inverse  is 


(I  -  P) 


-l  = 


3.  676  2.206 


5.  882  2.034 


.257 

.685 


7.519  2.532 

0  10.204 


Using  Eguation  2.1 


n  =  (35)  (1,  0,  0,  0) 


(3.676  2.206  .763 

0  5.882  2.034 

0  0  7.519 

0  0  0 


=  (128.66,  77.21,  26.705,  8.995)  . 

After  rounding  to  the  nearest  integer  the  S.S.S.V.  is: 
n  =  (129,  77,  27,  9) 

The  computer  printout  of  the  terminal  session  for 
this  example,  showing  the  stocks  for  years  10,  50  and  100, 

can  be  found  in  Appendix  A,  confirming  the  steady  state 
stock  vector  above-  As  the  printout  shows,  steady  state  is 
achieved  by  t  =  10. 


Example  <b) 

This  is  a  non-hierarchical,  three-grade  system. 

The  inital  stock  vector  is: 

n.  (0)  =  (300,  200,  100)  . 

The  transition  matrix  is: 

/ *  6 0  -  15  .05  \ 

P  =  j  .75  .20  0  1 

\.  05  0  .90  /  . 

The  total  recruitment  is: 

R(t)  =  100  ,  fixed  for  all  t>1  . 


The  recruitment  proportion  vector  is; 


Then 


I  -  P  = 


and  the  inverse  is 


(I  -  P)~*  = 


.4  -.15  -.05 

-.75  .8  0 

^-.  05  0  .10 


4.267  .8  2.  133 

4  2  2 

2.  133  .4  11.066 


Using  Equation  2.1 

4.  267 

n_  =  (100)  (.70,  .07,  .23)  f  4 


8 

2 


2.133 

2 


2.  133  .4  11.066 

=  (375.7,  79.2,  417.8)  - 

After  rounding  to  the  nearest  integer  the  S.S.S.V.  is; 


JL  =  (376,  79,  418)  . 


The  computer  printout  of  the  terminal  session  for 
this  example,  showing  the  stocks  for  years  10,  50  and  100, 
can  be  found  in  Appendix  A,  confirming  the  steady  state 
stock  vector  above.  As  the  printout  shows,  steady  state  is 
achieved  by  t  =  50. 


2.  Option  2 •  Additive  Recruitment 

This  option  addresses  the  fixed  additive  increase 
(decrease)  in  recruitment  for  each  time  period.  Here, 
R  ♦  (t- 1)  M,  replaces  R  (t)  in  Equation  1.2  resulting  in 


jn  (t)  =  n(t-  1)  P  ♦  Rr  ♦  (t-1)  Mr.  . 


(eqn  2.2) 


is  obtained.  The  steady  state  distribution  vector  is 
achieved  by  normalizing  ft,  which  results  in 

t  £  [I-*]"1  (eqn  2.3) 

tV  al' 

After  steady  state  has  been  reached  it  follows  that 
the  formula 

n(t)  =  tfl£  for  "large"  t, 

can  be  used  to  approximate  the  stocks  at  time  t.  The  two 
examples  from  section  Al  are  worked  out  below  to  illustrate 
Equation  2.3  . 

Example  (a) 

This  is  a  hierarchical,  four  category  system. 

The  inital  stccx  vector  is: 

_n(Q)  =  (129,  74,  28,  1  1) 

The  transition  matrix  is: 


The  total  recruitment  is: 

R  (t)  =  P.  «•  <t— 1 ) M  where  R  =  35  and  M  =  5  . 
The  recruitment  proportion  vector  is: 

r  =  (1,  0,  0,  0) 
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=  (.53,  .  32,  .  11  ,  .04) 

is  the  steady  state  distribution  vector.- 

The  computer  printout  of  the  terminal  session  for 
this  example,  showing  the  stocks  for  years  10,  50,  100, 

900,  and  975  can  be  found  in  Appendix  B,  confirming  that  in 
steady  state  there  will  be  53%,  32%,  11%,  and  4%  respec¬ 

tively,  in  categories  1,2,3,  and  4.  As  the  printout  shows, 
steady  state  distribution  (to  two  decimal  accuracy)  is 
achieved  oy  t  =  975  . 

Example  (b) 

This  is  a  non-hier archical,  three-grade  system. 

The  inital  stock  vector  is: 

n (0)  =  (300,  200,  100) 
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mj  rfrmrrp  m  vw 


The  transition  matrix  is: 


rw  wwwzrwsrx n  pgrwmoTny  *?i  *j\  vy\  v  k  TF5T  ^ITV,?71 


The  total  recruitment  is: 


R(t)  =  R  +  (t-1)M  where  R  =  100  and  M  =  25 
The  recruitment  proportion  vector  is: 

r  =  {.70,  .07,  .23)  . 


Then 


(-4  -.15 

-.75  .8 

-.05  0 

and  the  inverse  is 

(4.267  .8 

4  2 

2.133  .4 


-.05 

0 

.10 


2.  133 
2 

11.066 


Using  Equation  2.3 


n 


(.7,  .07,  .23) 
8.728 


✓  4.  267 
I  4.0 
\  2. 133 


.8 

2.0 


.4  11.066 


=  (  .  43,  .09,  .48) 

is  the  steady  state  distribution  vector. 

The  computer  printout  of  the  terminal  session 
this  example,  showing  the  stocks  for  years  10,  50, 


mr  nn  \ 


for 

100, 


25 


350,  and  400  can  be  found  in  Appendix  B,  confirming  that  in 
steady  state  there  will  be  43S,  9*,  and  48S  respectively,  in 
categories  1,2,  and  3.  As  the  printout  shows,  steady  state 
distribution  (to  two  decimal  accuracy)  is  achieved  by 
t  =  400  . 


3.  Option  3:  Multiplicative  fiecr uitment 

This  cption  addresses  the  multiplicative  increase 
(decrease)  in  recruitment  for  each  time  period.  Here, 
R0t_1  ,  replaces  R(t)  in  Eguation  1.2  resulting  in 

n(t)  =  n(t-1)P  ♦  Ee^-ir  .  (egn  2.4) 

The  values  that  6  can  assume  play  an  important  part  in  the 
derivation  of  the  steady  state  eguation.  First,  6  >0  must 
be  true,  otherwise  in  even  numbered  years  there  would  be 
negative  recruitment  (  R(t)<0  )  which  does  not  make  sense  in 
manpower  models.  There  are  three  cases  to  be  examined. 

Case  1  :  0<  0  <1 

Then 

limit  R (t)  =  0  , 

t  — ►  OO 

because  et_  1 — ►O  when  0<  €  <1  .  Therefore,  letting 

t— • ■  in  Equation  1.2  results  in 

n  =  nP  , 

or 

n  (I  -  P)  =  0-  . 

Since  1-p  has  a  unique  inverse  matrix  (see  Theorem  1)  ,  the 
only  solution  to  the  above  equation  is  n  =  0.  This  means 

that  such  a  system  becomes  empty  in  steady  state. 


Case  2: 


0  =  1 


Then 

limit  R  (t)  =  R  , 

t — ►  OO 

which  is  identical  to  option  1. 

Case  3:  0  >  1 

Assume  that, 

.  nm  - 

limit  - —  =  n 

t  —  OO  0*  “ 

exists.  This  amounts  to  assuming  that  for  "large"  t,  n(t) 
behaves  like  0*  £  .  Using  the  above  assumption,  Squation 
2.4  can  be  rewritten  as 

n  (t)  H(t-1)  R0*-* 

-  =  — - - -  p  +  - - -  r 

er  e*-»  e  9 


Taking  the  limit  as  t-*-<»  results  in 

^  ^  1  n  R 

n  =  n*  -  *P  ♦  -  r 

~  9  0  “, 


or 

n;  [ 01  -  P ]  =  Rr  . 

Finally,  knowing  that  the  matrix  0I-P  has  a  unique  inverse, 
since  0  >  1  (see  Theorem  1)  ,  the  result 

n  »  fir  [  9I-P  1~» 

is  obtained.  When  n  is  normalized,  the  resulting  vector  is 
n,  the  steady  state  distribution  vector  where. 


After  steady  state  has  been  reached  it  follows  that 
the  formula 


n(t)  =  9Tn.  for  “large"  t. 


can.  be  used  to  approximate  the  stocks  at  time  t.  Two  exam¬ 
ples  are  worked  out  below  to  illustrate  Equation  2.5 


Example  (a) 

This  is  a  hierarchical,  four  category  system. 


The  inital  stccx  vector  is: 


n<0)  =  (129,  74,  28,  1  1) 
The  transition  matrix  is: 


The  total  recruitment  is: 


h(t)  =  R6f-i  where  R  =  35  ,  0  =  1.01  . 

The  recruitment  proportion  vector  is: 

t  —  (1,0,  0,0)  • 


Then 


and  the  inverse  is 


[  01  -  P]-»  = 


546 

2.009 

.646 

.198 

0 

5.  556 

1.787 

.546 

0 

0 

6.993 

2.  137 

0 

0 

0 

9.259 

Using  Equation  2.5 


n  =  ^77 

nl ' 


(35)  (1,  0,  0,  0) 
223.965 


/  3.  546 

2.  009 

.646 

.198  \ 

/  0 

5.  556 

1.787 

.546  1 

0 

0 

6.99  3 

2.137  1 

\  o 

0 

0 

9.259  / 

=  (-  554,  .314,  .  101,  .031) 

is  the  steady  state  distribution  vector. 

The  computer  printout  of  the  terminal  session  for 
this  example,  showing  the  stocks  for  years  10,  50  ar.i  100, 

can  be  found  in  Appendix  C,  confirming  that  in  steady  state 
there  will  be  55*,  31*,  10%,  and  3%  respectively,  in  catego¬ 
ries  1,2,3,  and  4.  As  tae  printout  shows,  steady  state 
distribution  (to  two  decimal  accuracy)  is  achieved  by 
t  =  50  . 


Example  (b) 

This  is  a  non-hierarchical,  three-grade  system. 

The  inital  stock  vector  is: 

£(0)  =  (300,  200,  100) 

The  transition  matrix  is: 

(.60  -  15  .05  \ 

.75  .  20  0  I 

.05  0  .90  /  . 


The  total  recruitment  is: 


R(t)  =  Re*-*  where  R  =  100  ,  6  =  1.01 

The  recruitment  proportion  vector  is: 

r  =  (.70,  .07,  .23)  . 

Then 


ei  -  P  = 


and  the  inverse  is 

Tel  -  P]-* 

Using  Equation  2.5 


:4i 

-.15 

ID 

O 

• 

1 

.75 

.81 

0 

.05 

0 

.11 

03 

.75 

1.  33 

73 

1.92 

1.69 

83 

.  34 

9.92 

n  =  -t— 

-  n1» 


792.  15 


4.03 

-75 

1.83 

3.73 

1.92 

1.69 

1.83 

.34 

9.92 

=  (.44,  .09,  .47) 

is  the  steady  state  distribution  vector. 

The  computer  printout  of  the  terminal  session  for 
this  example,  showing  the  stocks  for  years  10,  50  and  100, 

can  be  found  in  Appendix  C,  confirming  that  in  steady  state 
there  will  be  44X,  9%,  and  47*  respectively,  in  categories 
1,2,  and  3.  As  the  printout  shows,  steady  state  distribu¬ 
tion  (to  two  decimal  accuracy)  is  achieved  by  t  =  100  . 


B.  TOTAL  STSTEH  SIZE  OPTIONS 


The  options  in  this  section  are  variations  of  Equation 
1.6  where  «  (t) ,  representing  the  change  in  total  system 
size,  is  the  term  that  will  take  on  different  forms.  As  a 
result,  the  steady  state  equations  for  these  options  are 
also  obtained  in  somewhat  different  forms.  In  the  next  two 
subsections  the  steady  state  equations  will  be  derived  for 
the  system  size  options. 

1.  Potion  4:  Additive  System  size 

This  option  addresses  the  additive  increase 
(decrease)  in  total  system  size  for  each  time  period.  In 
this  case,  M  (t)  =  B,  a  fixed  amount  of  increase  (decrease) 
per  time  period.  Rewriting  Equation  1.6  results  in 

n(t)  =  _n  (t-  1)  Q  ♦  ttr_  .  (eqn  2.6) 

There  are  3  cases  to  be  considered. 

Case  1:  B  <  0 

In  this  case  jobs  are  being  eliminated  by  a  fixed 
amount,  M,  each  year  which  implies  that  eventually  N(t)  =  0 
is  reached  in  a  finite  amount  of  time.  This  implies  that 
n  =  0_  is  the  steady  state  stock  vector.  This  means  that 
sucn  a  system  becomes  empty  in  steady  state. 

Case  2:  »  *  0 

Khen  B ( t)  =  0,  N  (t)  is  fixed  and  can  be  denoted  V. 
It  follows  that  Equation  1.6  can  be  written  as 


n  (t)  =  n'(t-1)Q  . 


Since  the 


limit 
t  —  » 


n(t) 


n 


(eqn  2.7) 


31 


exists  [Ref.  4:  pp.  40-43,  48-50  ],  Equation  2.7  can  be 

written  as 

n  =  n  Q  , 


n_  (I  -  Q)  =  0  . 

Because  the  first  column  of  i-q  is  a  linear  combination  of 
all  the  other  columns  the  inverse  of  I-Q  does  not  exist. 
So,  to  solve  for  n,  it  is  necessary  to  introduce  the  addi¬ 
tional  constraint  of  n1_'  =  N  [Ref.  5:  p.  80  ].  Replacing 
the  first  column  of  the  I-Q  matrix  with  a  column  of  ones, 
denoted  by  (I-Q)*,  and  replacing  the  first  element  of  0  with 
N,  denoted  0[*,  results  in  the  following  equation 

_n  (I-Q)  *  =  0*  . 

Taking  the  inverse  results  in 


n  =  0  *[  (I-Q)*  ]-» 


(egn  2.8) 


the  steady  state  stock  vector.  After  steady  state  is 
reached,  then  ja(t)  =  n,  for  "large"  t.  Two  examples  are 
worked  cut  below  to  illustrate  Equation  2-8  . 

Example  (a) 

This  is  a  hierarchical,  four  category  system. 

The  inital  stocx  vector  is: 

n_(0)  =  (129,  74,  28  ,  1  1) 


The  transition  matrix  is. 


$ 


The  0*  vector  is: 

0*  ~  (  242,  0,  0,  0  ) 
The  recruitment  proportion  vector  is: 

£  *  (1.  0,  0,  0)  . 

Then 


I  -  Q  = 


Replacing  the  first  column  vith  ones  gives 


(I  -  Q)*  = 


/' 

102 

0 

°  \ 

.  17 

-.046 

°  1 

1 ' 

0 

.  13  3 

-.033  J 

\1 

0 

0 

.  098  / 

and  the  inverse  is 


((I  -  Q)*)~l  = 


r  .53  3 

.320 

.111 

.037 

-4.582 

3.133 

1.084 

.  365 

-5.354 

-3.212 

-  6.408 

2.158 

-5.435 

-3.261 

-1.128 

9.  824 

Using  Equation  2.8 


=  (24  2,  0 


I  -4. 
'  °'  0)  (  -5. 


r  .533 

.320 

.111 

.037 

-4.582 

3.133 

1.084 

.365 

-5.354 

-3.212 

6.408  . 

2.  158 

-5.435 

-3.261 

-1.128 

9.824 

=  (128.99  ,  77.44  ,  26  .  86  ,  8.95)  . 

After  rounding  to  the  nearest  integer  the  S.S.S.V.  is: 
n  =  (129,  77,  27,  9)  . 

The  computer  printout  of  the  terminal  session  for 
this  example,  showing  the  stocks  for  years  10,  50  and  100, 

can  be  found  in  Appendix  D,  confirming  the  steady  state 
stock  vector  above.  As  the  printout  shows,  steady  state  is 
achieved  by  t  =  50 

Example  (b) 

This  is  a  non-hierarchical,  three-grade  system. 

The  inital  stock  vector  is: 

n  (0)  =  (300,  200,  100)  . 

The  transition  matrix  is; 


The  0*  vector  is: 

0* 


60 

.15  .  05  ) 

i 

75 

.20  0 

05 

0  .90  ) 

r 

» 

.  74 

.  1 64 

.096 

.  785 

.  2035 

.0115 

.  085 

.0035 

.9115 

00, 

0,  0  )  . 

The  recruitment  proportion  vector  is: 

r  =  (.70,  .07,  .23)  . 

Then 


I  -  Q  = 


26 

-.  164 

-.096 

785 

.  797 

-.012 

085 

-.  004 

.  089 

Replacing  the  first  column  with  ones  gives 


n 

-.  164 

-.096 

<1  -  Q)  *  =  (  1 

.797 

-.012 

\l 

-.004 

.089 

and  the  inverse  is 


((I  -  Q)*)-»  « 


.430  .090  .479  \ 

-.611  1.127  -.516  ) 

-4.888  -.981  5.869  ' 


Using  Equation  2.8 


n  =  (  6  00,  0,  0) 


.430  .090  .479 
-.611  1.127  -.516 
-4.888  -.981  5.869 


=  (  258.0,  54.4,  287.4)  . 

After  rounding  to  the  nearest  integer  the  S.S.S.V 
n  =  (258,  54,  287)  . 


The  computer  printout  of  the  terminal  session  for 
this  example,  shoving  the  stocks  for  years  10,  50  and  100, 

can  be  found  in  Appendix  Q,  confirming  tne  steady  state 
stock  vector  aoove.  As  the  printout  shovs,  steady  state  is 
achieved  by  t  =  50  . 

Case  3:  M  >  0 

Equation  1.6  can  then  be  rewritten  as 

n  (t)  n  (t- 1)  « 

N(t)  ~  N(t)  °  *  N  (t)  £  . 

Iterating  the  equation,  N(t)  =  N(t-1)  ♦  M,  for  t  =  1,2,... 
will  shew  that,  N  (t)  =  N  (0)  +  tM  for  all  t  =  1,2,...  and 

the  ref  ore, 

a  n ( t- 1  > 

limit  -  =0  and  limit  -  =  1 

t  — -  n (t)  t— •»  N(t-i)  ♦  a 

Then  writing 

n.  (t)  n (t-1 )  N  (t-1) _  3 

N  (t)  ~  N  (t-1 )  *[N(t-1)  +  d]*Q  +  N(t)£  , 

and  taking  the  limit  as  t— ►  oo  yields 

n  s  j  Q  ,  or  n_  (I-Q)  =  0_ 

where 

-  ..  Jl<t) 

n  =  limit  - 

t  —  oo  N  <t> 

is  the  steady  state  distribution  vector. 

This  is  identical  to  the  3  =  0  case,  except  here  the 
constraint  n1_*  =  1  must  hold.  Therefore,  by  replacing  the 
first  column  of  the  I-Q  matrix  with  a  column  or  ones. 


■J. . 


'J  's'.' 
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denoted  by  (I-Q)*,  and  replacing  the  first  element  of  0  with 
1,  denoted  J)**,  the  solution  is 

n_  =  0.**C  <I-Q)  *]~l  #  (eqn  2.9) 

the  steady  state  distribution  vector. 

After  steady  state  has  been  achieved,  Ji_(t)  can  be 
found  using 

n(t)  =  Qf  H  (0)  ♦  tS]n 

for  "large”  t.  Two  examples  are  worked  out  below  to  illus¬ 
trate  Equation  2.9  . 


Example  (a) 


This  is  a  hierarchical,  four  category  system. 


The  inital  stock  vector  is: 

n(0)  =  (129,  74,  28;  11) 
The  transition  matrix  is: 


and 


The  0**  vector  is: 


The  recruitment  proportion  vector  is 


r  =  (1*0,  0,0)  - 


Then 


Replacing  the  first  column  with  ones,  gives 


(I  -  Q)*  * 


102 

17 

0 

0 


and  the  inverse  is 


((I  -  Q)*)~l  = 


.533 

.320 

.111 

.  037 

-4.582 

3.  133 

1.084 

.  365 

-5.354 

-3.212 

6.408 

2.158 

-5.435 

-3.261 

-1.128 

9.  824 

Using  Eguation  2.9 


n  =  (1,  0,  0,  0) 


.533 

.320 

.111 

.  037 

-4.582 

3.  133 

1.084 

.  365 

-5.354 

-3.212 

6.408 

2.  158 

-5.435 

-3.261 

-1.128 

9.  824 

=  (.53*  .  32,  .  11  ,  .04) 


is  the  steady  state  distribution  vector. 

The  computer  printout  of  the  terminal  session  for 
this  example,  showing  the  stocks  for  years  10,  50,  100, 
900,  and  950  can  be  found  in  Appendix  D,  confirming  that  in 


steady  state  there  will  be  53%,  32%,  11%,  and  4%  respec 

tively,  in  categories  1,2,3,  and  4.  As  the  printout  shows 
steady  state  distribution  {to  two  decimal  accuracy)  i 
achieved  by  t  =  950 

Example  (b) 

This  is  a  non-hierarchical,  three-grade  system. 

The  inital  stock  vector  is: 


n_(0)  =  (300,  200,  100) 
The  transition  matrix  is: 


and 


i 

{  .60 

.  15 

.05  ^ 

F  = 

.75  . 

.  20 

\ 

^  .05 

0 

.90  / 

9 

/  *74 

m 

164 

.096 

Q  =  | 

.785 

m 

2035 

.0115 

\  .085 

m 

0035 

.9115 

The  0**  vector  is: 

0**  =  (  1,  0,  0,  )  . 
The  recruitment  proportion  vector  is: 

r  =  (.70,  -07,  .  23)  . 

Then 


26 

-.  164 

-.096  \ 

785 

.  797 

-.012 

085 

-.  004 

.089  / 

• 

Replacing  the  rirst  column  with  ones,  gives 


and  the  inverse  is 


( (I 


(.430 
-.611 
-4.888 


.091  .479 

1.127  -.516 

-.981  5.869 


Using  Equation  2.9 


(.430 
-.611 
-4.888 


.091  .479 

1.127  -.516 

-.  981  5.  869 


=  (  .  43.  .09,  .48) 

is  the  steady  state  distribution  vector. 

The  computer  printout  of  the  terminal  session  for 
this  example,  showing  the  stocks  for  years  10,  50,  100, 

350,  and  400  can  be  found  in  Appendix  D,  confirming  that  in 
steady  state  there  will  be  4  35?,  955,  and  48%  respectively,  in 
categories  1,2,  and  3.  As  the  printout  shows,  steady  state 
distribution  (to  two  decimal  accuracy)  is  achieved  by 
t  =  400  . 


2.  Option  5:  Multiplicative  System  Size 

This  section  addresses  the  option  of  multiplicative 
increase  (decrease)  in  total  system  size  for  eacn  time 
period.  Starting  with 

N  (t)  =  QN(t-l),  for  t  =  1,2,... 
it  follows  that 

N  (t)  =  9*^(0),  for  t  =  1,2,... 
from  which  M  (t)  can  be  determined  as 


M  (t)  =  N  (t)  -  N  (t-1) 


(9-  1)  0r-»N  (0)  . 


Then,  letting  <0- 1)  S**-1  N  (0)  replace  M (t)  in  Equation  1.6, 
results  in 

a(t)  =  &(t-1)Q  ♦  [  (0-1 )  0^” 1  N  (0 )  ]  £  .  (eqn  2.  10) 

First,  6  >  0  must  be  true  otherwise  negative  or  zero  system 
size  may  result.  This  leaves  two  cases  to  be  examined. 

Case  Is  0  <  6  <  1 

In  this  case 

limit  H  ( t)  =0 
t— ®o 

because  as  t  — <*>  ,  e*-i— ►  o  for  0  <  0  <1  .  Finally,  if 
9  =  1  then  e-1  =  0  and  N  (t)  =  0  for  all  t.  Therefore, 

n_  =  limit  n(t)  =  0 
t  -*■  «• 

must  also  hold.  This  means  that  such  a  system  becomes  empty 
in  steady  state. 

Case  2:  0  >  1 

Equation  2. 10  can  be  rewritten  as 

n  (t)  _  a  (t-1 )  (9-1)9*-»H(0) 

N(t)  ~  N(t)  Q  +  N(t)  “  . 

Since  N  (t)  can  be  written  as  0N(t-1)  or  as  6rN(0),  the  above 
equation  gives 

n  (t)  n  ( t-  1 )  1  (0-1)  9*-» N  (0) 

-  =  - •*  -  «q  *■  - 1 -  r 

N  (t)  N  (t- 1  )  0  0*N(O) 

or 

n  (t)  n  ( t- 1 )  1  0-1 

— - -  =  - •  -  *Q  ♦  - r 

N  (t)  N  (t- 1  )  0  0  “ 

— <*>  results  in 


Taking  the  limit  as  t 
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n 


limit 
t  -*•  •« 


n  (t) 

H  (t) 


is  the  steady  state  distribution  vector. 

Since  e  >  1  in  this  case,  the  matrix  ei-Q  has  a 
unique  inverse  according  to  Theorem  2,  so  the  result 

n  =  (©-1)  r[  ei-Q]-1  (eqn  2.11) 

is  obtained  for  the  steady  state  distribution  vector.  After 
steady  state  has  been  reached  it  follows  that  the  formula 

n_(t)  =  6rN  (0)  n 

can  be  used  tc  find  the  stocks  at  time  t.  Two  examples  are 
worked  out  below  to  illustrate  Equation  2-11  . 

Example  (a) 

This  is  a  hierarchical,  four  category  system. 

The  inital  stcck  vector  is: 

n (0)  =  (129,  74,  28  ,  1  1) 


The  transition  matrix  is: 


(] 


0 

0 

.033 

.902 


The  multiplicative  factor  is: 

9  =  1.01  . 

The  recruitment  proportion  vector  is: 

r  =  (1,  0,  0 ,  0)  . 

Then 

(.112  -. 102 

-.124  .18 

1  0 

-.098  0 

and  the  inverse  is 


(91  -  Q) 


0 

0 

-.033 

.  108 


55.4  1 

31.40 

10.  10 

3.09 

51.04 

34.48 

11.09 

3.39 

50.35 

28.53 

16.  17 

4.94 

50.  28 

28.49 

9.  17 

12.06 

Using  Equation  2.11 


n  =  (.0  1)  (1 ,0,0,0) 


=  (.  554,  .  3  14,  .  101,  .031) 
is  the  steady  state  distribution  vector. 


55.4  1 

31.40 

10.  10 

3.09 

51.04 

34.48 

11.  09 

3.  39 

50.35 

28.53 

16.  17 

4.94 

50.  28 

28.49 

9.  17 

12.06 

The  computer  printout  of  the  terminal  session  for 
this  example,  shoving  the  stocks  for  years  10,  50,  and  100, 
can  be  found  in  Appendix  E,  confirming  that  in  steady  state 
there  will  be  55%,  31%,  10%,  and  3%  respectively,  in  catego¬ 
ries  1,2,3,  and  4.  As  the  printout  shows,  steady  state 
distribution  (to  two  decimal  accuracy)  is  achieved  by 
t  =  50  . 

ixaorgle  (b) 

This  is  a  non-hierarchical,  three-grade  system. 

The  inital  stock  vector  is: 

n  (0)  *  (300,  200,  100)  . 

The  transitior  matrix  is: 

/  * 60  -  15  .05  \ 

P-  *75  *2°  °  ) 

\  • 05  0  .90  /  , 

and 


/-7“ 

a  *  I  -785 
V  .  085 


.164  .096  \ 

.2035  .0115  ] 

.0035  .9115  / 


The  multiplicative  factor  is: 


6  =  1.01  . 


The  recruitment  proportion  vector  is: 


and  the  inverse  is 


(61  -  Q)-» 


(  45.  32  9.41  45.  27 

44.69  10.52  44.79 

40.  70  8.50  50.  81 


Using  Equation  2.11 


n  =  (.0  1)  (.7, .07, .23) 


9.41  45.27  \ 

10.52  44.79  J 

8.  50  50.81  / 


=  (.44,  .09,  .47) 

is  the  steady  state  distribution  vector. 

The  computer  printout  of  the  terminal  session  for 
this  example,  showing  the  stocks  for  years  10,  50,  and  100, 
can  be  found  in  Appendix  E,  confirming  that  in  steady  state 
there  will  oe  44?,  9 %,  and  47?  respectively,  ir.  categories 
1,2,  and  3.  As  the  printout  shows,  steady  state  distribu¬ 
tion  (tc  two  decimal  accuracy)  is  achieved  by  t  =  100  . 


This  chapter  is  a  synopsis  of  the  options  that  describe 
the  five  different  conditions  of  growth.  In  all  options, 
the  restriction  of  wt  >  0  for  all  i  =  1,2,...,k  was  assumed, 
since  attrition  must  be  allowed  from  any  category  in  a 
manpower  system.  Options  1  through  3  employed  the  basic 
prediction  equation  (Equation  1.2)  and  used  constant,  addi- 
tively  and  multiplicatively  increasing  (decreasing)  recruit¬ 
ment,  R(t),  to  derive  the  steady  state  behavior  of  a 
manpower  system.  For  options  4  and  5,  Equation  1.6  used 
additively  and  multiplicatively  increasing  (decreasing) 
total  system  size,  M(t),  to  derive  the  steady  state  behavior 
of  a  manpower  system. 

Derivation  of  the  steady  state  stock  vector  or  distribu¬ 
tion  vector  required  some  standard  results  of  matrix  algebra 
summarized  in  Theorems  1  and  2  of  Chapter  II.  In  addition 
to  the  analytical  derivation  of  the  steady  state  equations, 
examples  have  been  worked  out  in  each  case  (excepting  some 
trivial  cases)  to  illustrate  the  use  of  the  analytical 
results.  The  formulas  have  been  programmed  in  the  AFL 
language  and  the  program  listing  is  given  in  Appendix  F. 
This  program  (called  a  function  in  APL)  has  been  integrated 
into  the  main  APL  program  at  NPS  which  forecasts  future 
stocks  of  a  manpower  system  using  the  Markov  chain  theory 
summarized  in  Chapter  I.  Appendices  A  through  E  show  the 
printouts  of  the  terminal  sessions  where  the  examples 
referred  to  above  are  worked  out  to  verify  the  steady  state 
results  obtained  in  each  analytical  example. 


APPENDIX  A 

COMPUTER  PBI1TOUTS  OF  EXAMPLES  FOB  OPTION  1 


1.  Example  (a) 


START 

DO  YOU  WISH  TC  ENTER  DATA? 

◦  NC 

1  YES 

* 

1 

ENTEfi  THE  NDMEEE  OF  THE  MODEL  TYPE 

1  MARKOV  HIERARCHICAL 

2  MARKOV  LENGTH  OF  SERVICE 

3  MARKOV  GENERAL 

4  VACANCY 

1 

ENTER  N  (INITIAL  STOCK  VECTOR) 

129  74  28  11 

ENTER  THE  PROMOTION  RATE  VECTOR. 

THIS  VECTOR  SHOULD  INCLUDE  THE  PROMOTION  RATES  FOR  THE  FIRST 
3  CLASSES. 

.  102  .046  .033 

ENTER  THE  HASTAGE  RATE  VECTOR. 

THE  4TH  VALUE  SHCULD  INCLUDE  THE  TOTAL  RATE  FOR  THAT 
STATE  DUE  TO  EITHER  PROMOTION  OR  HASTAGE. 

.17  .124  .1  .098 

ENTER  THE  NUMBER  OF  THE  RECRUIT  TYPE 

1  FIXED  RECRUIT  VECTOR 

2  ADDITIVE  (RECRUIT  SIZE) 

3  MULTIPLICATIVE  (RECRUIT  SIZE) 

4  ADDITIVE  (SYSTEM  SIZE) 

5  MULTIPLICATIVE  (SYSTEM  SIZE) 

1 

ENTER  R  (RECRUITMENT  VECTOR) 

35  0  0  0 

ENTER  THE  PERCENT  CODE 
0  NO  GRADE  PERCENTAGES 

1  GRADE  SIZE  AS  PERCENT  OF  TOTAL  GRADE  SIZE 

2  GRADE  SIZE  AS  PERCENT  OF  ORIGINAL  GRADE  SIZE 

7  QUIT  PROGRAM 

1 

WOULD  YOU  LIKE  TO  SEE  THE  ENTERED  DATA? 

0  NO 

1  YES 

1 


47 


0  0 

0.046  0 
0.867  0.033 
0  0.902 


WOULD  YOU  LIKE  TO  CHANGE  ANY  OF  THE  DATA? 

0  NO 
1  YES 

7  QUIT  PBC3EAM 
0 

ENTER  THE  NUMEER  OF  THE  YEAR  YOU  WISH  TO  SEE 
:  10 

DO  YOU  WISH  TO  SEE  THE  INTERVENING  YEARS? 

0  NO 
1  YES 

« 

0 


P  HA IRIX 
0.728  0.102 
0  0.83 

0  0 

0  0 

N  VECTOR 
129  74  28  1  1 
OPTION  =1 
R  VECTOR 
35  0  0  0 


T  N  PERCENT  R 


0  1 

129 

53) 

2 

74 

31) 

3 

28 

12) 

4 

11 

5) 

TOTAL 

242 

100) 

10  1 

2 

129 

77 

II 

3 

27 

11 

4 

10 

4 

TOTAL 

242 

100) 

DO  YOU  KISH  TO  SEE  ANY  OTHER  YEARS? 

0  NC 
1  Y  ES 

7  QUIT  PROGRAM 

« 

1 

ENTER  THE  NUMEER  OF  THE  YEAR  YOU  WISH  TO  SEE 
50 

DO  YOU  WISH  TC  SEE  THE  INTERVENING  YEARS? 

0  NO 
1  YES 

0 


1 

129 

[  53) 

2 

77 

32 

3 

27 

:  ii 

4 

9 

4' 

TOTAL 

242 

[  100 

35 


DO  YOU  WISH  TC  SEE  ANY  OTHER  YEARS? 

0  HO 
1  YES 

7  QUIT  PROGRAM 

1 

ENTER  THE  NUMEER  OF  THE  YEAR  YOU  WISH  TO  SEE 
100 

DO  YOU  WISH  TO  SEE  THE  INTERVENING  YEARS? 

0  NO 
1  YES 


100 


1 

129 

[  53) 

2 

77 

32 

3 

27 

11 

4 

9 

i  4 

TOTAL 

242 

[  100 

35 


DO  YOU  WISH  TO  SEE  ANY  OTHER  YEARS? 
0  NO 
YES 

QUIT  PROGRAM 


1 

7 


DO  YOU  WISH  TO  SEE  THE  STEADY  STATE  VECTOR? 
0  NO 
1  >  ES 


1 


999 


1 

129 

53) 

2 

77 

1  32 

3 

27 

11 

4 

9 

4 

TOTAL 

242 

100) 

35 


ARE  YOU  THROUGH? 
0  NO 
1  YES 


1 


2.  Exa  a  pie  (b) 


START 


DC  YOU  KISH  TC  ENIER  DATA? 
0  NC 
YES 


1 


1 


ENTER 

1 

2 

3 

4 


THE  NUMEER  OF  THE  MODEL  TYPE 
MARKOV  HIERARCHICAL 
MARKOV  LENGTH  OF  SERVICE 
MARKOV  GENERAL 
VACANCY 


ENTER  N  (INITIAL  STOCK  VECTOR) 

300  200  100 

ENTER  P  (TRANSITION  MATRIX)  BY  SONS 
ENTER  1TH  ROW 

.6  .15  .05 
ENTER  21H  ROW 

.75  .2  0 
ENTER  3TH  ROW 

.  05  0  .9 

ENTER  THE  NUMBER  OF  THE  RECRUIT  TYPE 

1  FIXED  RECRUIT  VECTOR 

2  ACDITIV1:  (RECRUIT  SIZE) 

3  MULTIPLICATIVE  (RECRUIT  SIZE) 

4  ADDITIVE  (SYSTEM  SIZE) 

5  MULTIPLICATIVE  (SYSTEM  SIZE) 

• 

1 

ENTER  R  (RECRUITMENT  VECTOR) 

70  7  23 

ENTER  THE  PERCENT  CODE 
0  NO  GRADE  PERCENTAGES 

1  GRADE  SIZE  AS  PERCENT  OF  TOTAL  GRADE  SIZE 

2  GRADE  SIZE  AS  PERCENT  OF  ORIGINAL  GRADE  SIZE 

7  QUIT  PROGRAM 

1 

WOULD  YOU  LIKE  TO  SEE  THE  ENTERED  DATA? 

0  NO 

1  YES 

1 


P  MATRIX 
0.6  0.15  0.05 

0.  75  0.2  0 

0.  05  0  O.S 
N  VECTOR 
300  200  100 
OPTION  =1 
R  VECTOR 
70  7  23 

WOULD  YOU  LIKE  TO  CHANGE  ANY  CF  THE  DATA? 
0  NO 
1  YES 

7  QUIT  PRCGRAM 


ENTER  THE  NUMBER  OF  THE  YEAR  YOU  WISH  TO  SEE 


DC  YOU  WISH  TC  SEE  THE  INTERVENING  YEARS? 
0  NO 
1  YES 


50 


.*  */  J*  ‘  *  W  “  •  ‘  •  ‘  *  "  «  “  -  '.*•*.**  •*»  .*  /'  »"•  •*'  .**  »  .  • 

•>  .>,  j'.v'.VVvV.-  "-■v'-'t’. lV.V 


T 

N 

PERCENT 

0 

1 

300  | 

1  50) 

2 

200  { 

33) 

3 

100  i 

17) 

TOTAL 

600  | 

[  100) 

10 

1 

351  | 

[  48 

2 

75  | 

10 

3 

302  | 

41 

TOTAL 

728  | 

121) 

100 


DO  YOU  NISH  TC  SEE  ANY  OTHER  YEARS? 

0  NO 
1  YES 

7  QUIT  PROGRAM 

1 

ENTER  THE  NUMEER  OF  THE  YEAR  YOU  PISH  TO  SEE 
50 

DO  YOU  WISH  TC  SEE  THE  INTERVENING  YEARS? 

0  NO 
1  YES 


50  1 

375  ( 

43 

2 

79 

9 

3 

415 

48 

TOTAL 

868 

145 

100 


DO  YOU  WISH  TC  SEE  ANY  OTHER  YEARS? 

0  NO 
1  YES 

7  QUIT  PROGRAM 

1 

ENTER  THE  NUMEER  OF  THE  YEAF  YOU  WISH  TO  SEE 
1  00 

DO  YOU  WISH  TC  SEE  THE  INTERVENING  YEARS? 

0  NO 
1  YES 


100 


100  1 

376  ( 

[  43 

2 

75  | 

9 

3 

418  | 

48 

TOTAL 

873  i 

146 

DC  YOU  WISH  TC  SEE  ANY  OTHER  YEARS? 

0  NO 
1  YES 

7  QUIT  PROGRAM 
0 

DO  YOU  WISH  TC  SEE  THE  STEADY  STATE  VECTOR? 
0  NO 
1  YES 


1 


51 


APPEBDIX  B 

COMPOTES  PRIHTOOTS  OF  EXAMPLES  FOB  OPTION  2 


1  *  Exa  mcle  (a) 


START 

DO  YOO  WISH  TO  ENTER  DATA? 

0  NC 

1  YES 

«* 

1 

ENTER  THE  NOMEER  OF  THE  MODEL  TYPE 

1  MARKOV  HIERARCHICAL 

2  MARKOV  LENGTH  OF  SERVICE 

3  MARKOV  GENERAL 

4  VACANCY 

1 

ENTER  N  (INITIAL  STOCK  VECTOR) 

129  74  28  11 

ENTER  THE  PROMOTION  RATE  VECTOR. 

THIS  VECTOR  SHOULD  INCLUDE  THE  PROMOTION  RATES  FOR  THE  FIRST 
3  CLASSES. 

.  102  .0  46  .033 

ENTER  THE  WASTAGE  RATE  VECTOR. 

THE  4TH  VALUE  SHOULD  INCLUDE  THE  TOTAL  RATE  FOR  THAT 
STATE  DUE  TO  EITHER  PROMOTION  OR  WASTAGE. 

.17  .124  .1  .098 

ENTER  THE  NUMBER  OF  THE  RECRUIT  TYPE 

1  FIXED  RECRUIT  VECTOR 

2  ADDITIVE  (RECRUIT  SIZE) 

3  MULTIPLICATIVE  (RECRUIT  SIZE) 

4  ADDITIVE  (SYSTEM  SIZE) 

5  MULTIPLICATIVE  (SYSTEM  SIZE) 

2 

ENTER  R  (RECRUITMENT  VECTOR) 

35  0  0  0 

ENTER  ADDITIVE  INCREASE 
* 

5 

ENTER  THE  PERCENT  CODE 
0  NO  GRADE  PERCENTAGES 

1  GRADE  SIZE  AS  PERCENT  OF  TOTAL  GRADE  SIZE 

2  GRADE  SIZE  AS  PERCENT  OF  ORIGINAL  GRADE  SIZE 
7  QUIT  PROGRAM 

1 


53 


WOULD  YOU  LIKE  TO  SEE  THE  ENTERED  DATA? 
0  NO 


P  MATRIX 

0.  728  0.102  0  0 

0  0.83  0.  046  0 

0  0  0.  867  0.033 

0  0  0  0.902 

N  VECTOR 
129  74  28  1  1 

OPTION  =2 
INC  =5 

TOTAL  RECRUITMENT 
35 

RECRUITMENT  PROPORTION  VECTOR 
10  0  0 

WOULD  YCU  LIKE  TO  CHANGE  ANY  OF  THE  DATA? 
0  NO 
1  YES 

7  QUIT  PROGRAM 


ENTER  THE  NUMEER  OF  THE  YEAR  YOU  WISH  TO  SEE 
1  0 

DO  YOU  WISH  TO  SEE  THE  INTERVENING  YEARS? 

0  NO 
1  YES 


PERCENT 


1 

2 

3 

4 

TOTAL 


1 

2 

3 

4 

TOTAL 


129 

531 

74 

31) 

28 

12 

1  1 

5 

242 

100) 

248 

63) 

106 

27) 

29 

7) 

10 

3 

392 

162 

DO  YOU  WISH  TC  SEE  ANY  OTHER  YEARS? 
0  NC 
1  Y  ES 

7  QUIT  PROGRAM 


ENTER  THE  NUMEER  OF. THE  YEAR  YOU  WISH  TO  SEE 
50 

DO  YCU  WISH  TO  SEE  THE  INTERVENING  YEARS? 

0  NO 


1 

980 

58 

2 

523 

31 

3 

152 

9 

4 

39 

2 

TOTAL 

1695 

700 

DO  YOU  WISH  TO  SEE  ANY  OTHER  YEARS? 
0  NO 
1  YES 

7  QUIT  PROGRAM 


ENTER  THE  NUMEER  OF  THE  YEAR  YOU  WISH  TO  SEE 
1  00 

DO  YOU  WISH  TO  SEE  THE  INTERVENING  YEARS? 

0  NO 
1  YES 


1 

2 

3 

4 

TOTAL 


1899 

1075 

343 

102 

3419 


56} 

31) 

10 

3) 

1413 


DO  YOU  WISH  TC  SEE  ANY  OTHER  YEARS? 
0  NC 
1  YES 

7  QUIT  PROGRAM 


ENTER  THE  NUMEER  OF  THE  YEAR  YOU  WISH  TO  SEE 
9  00 


DO  YOU  WISH  TO 

SEE  THE  INTERVENING 

YEARS? 

0  NO 

1  YES 

• 

0 

900  1 

16605 

54) 

2 

9898 

32) 

3 

3395 

11} 

4 

1130 

4) 

TOTAL 

31028 

12821) 

4530 

DC  YOU  WISH  TC  SEE  ANY  OTHER  YEARS? 
0  NC 
1  Y  ES 

7  QUIT  PROGRAM 


ENTER  THE  NUMEER  OF  THE  YEAR  YOU  WISH  TO  SEE 
9  75 

DO  YOU  WISH  IC  SEE  THE  INTERVENING  YEARS? 

0  NO 
1  YES 


0 


r-VV1 


1 

2 

3 

4 

TOTAL 


17984 

10725 

3681 

1226 

33617 


53’ 

32 

11( 

4j 

13891 


DO  YOU  WISH  TC  SEE  ANY  OTHER  YEARS? 
0  NO 
1  YES 

7  QUIT  PROGRAM 


DO  YOU  NISH  TC  SEE  THE  STEADY  STATE  VECTOR? 
0  NC 


PERCENTAGES  ARE  IN  STEADY  STATE 


1 

2 

3 

4 

TOTAL 


17984 

10725 

3681 

1226 

33617 


ARE  YOU  THROUGH? 
0  NO 
1  YES 


2-  Exa  mple  <b) 


START 


DO  YOU  NISH  TC  ENTER  DATA? 
0  NO 
1  YES 


ENTER  THE  NUHEEP  OF  THE  MODEL  TYPE 

1  MARKOV  HIERARCHICAL 

2  MARKOV  LENGTH  OF  SERVICE 

3  MARKOV  GENERAL 

4  VACANCY 


ENTEP  N  (INITIAL  STOCK  VECTOR) 

•  300  200  100 

ENTER  P  (TRANSITION  MATRIX)  BY  ROWS 
ENTER  1  TH  RON 

.6  .15  .05 

ENTER  2 TH  ROW 

.75  .2  0 

ENTER  3  TH  ROW 


ENTER  THE  NUMEER  OF  THE  RECRUIT  TYPE 

1  FIXED  R EC R 01 T  VECTOR 

2  ADDITIVE  (RECRUIT  SIZE) 

3  MULTIPLICATIVE  (RECRUIT  SIZE) 

4  ADDITIVE  (SYSTEM  SIZE) 

5  MULTIPLICATIVE  (SYSTEM  SIZE) 

2 

ENTER  fi  (RECRUITMENT  VECTOR) 

70  7  23 

ENTER  ADDITIVE  INCREASE 
25 

ENTER  THE  PERCENT  CODE 
0  NO  GRADE  PERCENTAGES 

1  GRADE  SIZE  AS  PERCENT  OF  TOTAL  GRADE  SIZE 

2  GRADE  SIZE  AS  PERCENT  OF  ORIGINAL  GRADE  SIZE 

7  QUIT  PROGEAM 

1 

WOULD  YOU  LIKE  TO  SEE  THE  ENTERED  DATA? 

0  NO 

1  YES 

1 


TOTAL  RECRUITMENT 

100 


RECRUITMENT  PROPORTION  VECTOR 
0.7  0.07  0.  23 

WOULD  YOU  LIKE  TO  CHANGE  ANY  OF  THE  DATA? 

0  NC 
1  YES 

7  QUIT  PROGRAM 
0 

ENTER  THE  NUMEER  OF  THE  YEAR  YOU  WISH  TO  SEE 
1  0 


57 


ENTER  THE  NUMEER  OF  THE  YEAR  YOD  NISH  TO  SEE 


DO  YOU  NISH  TO  SEE  THE  INTERVENING  YEARS? 
0  NO 
1  YES 


1 

2 

3 

TOTAL 


32550 

6853 

35585 

74987 


43 

9 

47 

12498 


8825 


DO  YOU  WISH  TC  SEE  ANY  OTHER  YEARS? 
0  NO 
1  YES 

7  QUIT  PROGRAM 


ENTER  THE  NUMEER  OF  THE  YEAR  YOU  FISH  TO  SEE 
4  00 


DO  YOU  NISH  TC  SEE  THE  INTERVENING  YEARS? 
0  NO 
1  YES 


1 

2 

3 

TOTAL 


37246 

7843 

40608 

85897 


%3 

46 

14316 


10075 


DO  YOU  NISH  TC  SEE  ANY  OTHER  YEARS? 
0  NO 
1  YES 

7  QUIT  PROGRAM 


DO  YOU  NISH  TC  SEE  THE  STEADY  STATE  VECTOR? 
0  NC 
1  Y  ES 


PERCENTAGES  ARE  IN  STEADY  STATE 


999  1 

37246 

43) 

2 

7843  r 

9) 

3 

40808  i 

48) 

TOTAL 

85897 

100) 

10075 

ARE  YOU  THROUGH? 
0  NC 
1  YES 


^  fL*r-  ^ 


APP ENDII  C 

COMPOTES  PRINTOUTS  OF  EXAMPLES  FOB  OPTION  3 


1-  Example  (a) 


START 

DO  YOU  WISH  TO  ENTER  DATA? 

0  NC 

1  YES 

« 

1 

ENTER  T BE  NUMEER  OF  THE  MODEL  TYPE 

1  MARKOV  HIERARCHICAL 

2  MARKOV  LENGTH  OF  SERVICE 

3  MARKOV  GENERAL 

4  VACANCY 

1 

ENTER  N  (INITIAL  STOCK  VECTOR) 

129  74  28  11 

ENTER  THE  PROMOTION  RATE  VECTOR. 

THIS  VECTOR  SHOULD  INCLUDE  THE  PROMOTION  RATES  FOR  THE  FIRST 
3  CLASSES. 

.  102  .0  46  .033 

ENTER  THE  WASTAGE  RATE  VECTOR. 

THE  4TH  VALUE  SHOULD  INCLUDE  THE  TOTAL  RATE  FOR  THAT 
STATE  DUE  TO  EITHER  PROMOTION  OR  WASTAGE. 

.17  .124  .1  .098 

ENTER  the  NUMEER  OF  THE  RECRUIT  TYPE 

1  FIXED  RECRUIT  VECTOR 

2  ADDITIVE  (RECRUIT  SIZE) 

3  MULTIPLICATIVE  (RECRUIT  SIZE) 

4  ADDITIVE  (SYSTEM  SIZE) 

5  MULTIPLICATIVE  (SYSTEM  SIZE) 

3 

ENTER  R  (RECRUITMENT  VECTOR) 

35  0  0  0 

ENTER  MULTIPLICATIVE  FACTOR 

1  .01 

ENTER  THE  PERCENT  CODE 
0  NO  GRADE  PERCENTAGES 

1  GRADE  SIZE  AS  PERCENT  OF  TOTAL  GRADE  SIZE 

2  GRADE  SIZE  AS  PERCENT  OF  ORIGINAL  GRADE  SIZE 
7  QUIT  PROGRAM 


WOULD  ICO  LIKE  TC  SEE  THE  ENTERED  DATA? 
0  NO 

1  YES 


1 


0  0 

0.046  0 
0.  867  0.033 
0  0.902 


P  MATRIX 
0.  728  0.102 
0  0.83 

0  0 
0  0 
N  VECTOR 
129  74  28  1  1 
OPTION  =3 
FACT  =1.0  1 

TOTAL  RECRUITMENT 
35 

RECRUITMENT  PROPORTION  VECTOR 
10  0  0 


WOULD  YOU  LIKE  TO  CHANGE  ANY  CF  THE  DATA? 

0  NO 
1  YES 

7  QUIT  PROGRAM 
0 

ENTER  THE  NUMEER  OF  THE  YEAR  YOU  WISH  TO  SEE 
1  0 

DO  YOU  WISH  TC  SEE  THE  INTERVENING  YEARS? 

0  NO 
1  YES 

0 


T  N  PERCENT 


0 

1 

2 

3 

4 

TOTAL 

129 

74 

28 

1 1 
242 

J? 

;  >f; 

100 

10 

1 

137 

541 

2 

79 

31 

3 

27 

11} 

4 

10  1 

4) 

TOTAL 

253  i 

105 

R 


38 


DO  YOU  WISH  TO  SEE  ANY  OTHER  YEARS? 

0  NC 
1  YES 

7  QUIT  PRCG RAM 

• 

1 

ENTER  THE  NUMEER  OF  THE  YEAR  YOU  WISH  TO  SEE 
50 

DO  YOU  WISH  TO  SEE  THE  INTERVENING  YEARS? 

0  NO 
1  YES 

m 

0 


ol 


1 

2 

3 

4 

TOTAL 


DO  YOU  WI  Sti  TO  SEE  ANY  OTHER  YEARS? 
0  NO 
1  YES 

7  QUIT  PROGRAM 


ENTER  THE  NUMBER  OF  THE  YEAR  YOU  WISH  TO  SEE 
1  00 

DO  YOU  WISH  TC  SEE  THE  INTERVENING  YEARS? 

0  NO 
1  YES 


1 

2 

3 

4 

TOTAL 


DO  YOU  WISH  TC  SEE  ANY  OTHER  YEARS? 
0  NC 
1  YES 

7  QUIT  PROGRAM 


DO  YOU  WISH  TO  SEE  THE  STEADY  STATE  VECTOR? 
0  NC 
1  Y  ES 


PERCENTAGES  ARE  IN  STEADY  STATE 
999  1  336 

2  190 

3  61  i 

4  19  r 

TOTAL  606  T  1 


ARE  YOU  THROUGH? 
0  NO 
1  YES 


2.  Example  (b) 

SIART 

DO  YOU  KISH  TO  ENTER  DATA? 


CHU1 


ENTEfi  THE  NUMEEB  OF  THE  MODEL  TYPE 

1  MARKOV  HIERARCHICAL 

2  MARKOV  LENGTH  OF  SERVICE 

3  MARKOV  GENERAL 

4  VACANCY 

3 

ENTER  N  (INITIAL  STOCK  VECTOR) 

300  200  100 

ENTER  P  (TRANSITION  MATRIX)  BY  ROWS 
ENTER  1  TH  ROW 

.6  .15  .05 
ENTER  2TH  ROM 

.75  .2  0 
ENTER  3 TH  ROW 

.  05  0  .9 

ENTER  THE  NUMBER  OF  THE  RECRUIT  TYPE 

1  FIXED  RECRUIT  VECTOR 

2  ADDITIVE  (RECRUIT  SITE) 

3  MULTIPLICATIVE  (RECRUIT  SIZE) 

4  ADDITIVE  (SYSTEM  SIZE) 

5  MULTIPLICATIVE  (SYSTEM  SIZE) 

• 

3 

ENTER  R  (RECRUITMENT  VECTOR) 

70  7  23 

ENTER  MULTIPLICATIVE  FACTOR 

1.01 

ENTER  THE  PERCENT  CODE 
0  NO  GRADE  PERCENTAGES 

1  GRADE  SIZE  AS  PERCENT  OF  TOTAL  GRADE  SIZE 

2  GRADE  SIZE  AS  PERCENT  OF  ORIGINAL  GRADE  SIZE 

7  QUIT  PROGRAM 

1 

WOULD  YOU  LIKE  TO  SEE  THE  ENTERED  DATA? 

0  NO 

1  YES 

1 


0.6  0.15  0.  05 

0.  75  0.2  0 

0.05  0  0.9 

N  VECTOR 
300  200  100 

OPTION  =3 

FACT  =1.01 

TOTAL  RECRUITMENT 

100 

RECRUITMENT  PROPORTION  VECTOR 
0.7  0.07  0.  23 


63 


WOULD  YOU  LIKE  TO  CHANGE  ANY  OF  THE  DATA? 

0  NO 
1  YES 

7  QUIT  PHOGBAM 

* 

0 

ENTER  THE  NUMEER  OF  THE  YEAR  YOU  WISH  TO'  SEE 
1  0 

DO  YOU  WISH  TO  SEE  THE  INTERVENING  YEARS? 

0  NO 
1  YES 

0 

T  N  PERCENT  R 


0 

1 

300  | 

50 

2 

200  | 

;  33 

3 

100  l 

17 

TOTAL 

600  \ 

100) 

10 

1 

49) 

2 

78  1 

10 

3 

312  l 

41 

TOTAL 

761  | 

127 

DO  YOU  WISH  TO  SEE  ANY  OTHER  YEARS? 

0  NC 
1  YES 

7  QUIT  PROGRAM 

1 

ENTER  THE  NUMEER  OF  THEf  YEAR  YOU  WISH  TO  SEE 
50 

DO  YOU  WISH  TO  SEE  THE  INTERVENING  YEARS? 

0  NO 
1  YES 

m 

0 

50  1 

2 
3 

TOTAL 


DO  YOU  WISH  TC  SEE  ANY  OTHER  YEARS? 

0  NC 
1  YES 

7  QUIT  PROGRAM 

1 

ENTER  THE  NUMEER  OF  THE  YEAR  YOU  WISH  TO  SEE 
1  00 


575  | 

44 

121  | 

9 

603  | 

46 

1298  { 

216 

DO  YOU  WISH  TC  SEE  THE  INTERVENING  YEARS? 


APPENDIX  D 

COHPUTEB  PRINTOUTS  OF  EXAMPLES  FOB  OPTION  4 


1.  Case  J  Example  (a) 


START 

DO  YOU  WISH  TC  ENTER  DATA? 

0  NO 

1  YES 

1 

ENTER  THE  NUMBER  OF  THE  MODEL  TYPE 

1  MARKOV  HIERARCHICAL 

2  MARKOV  LENGTH  OF  SERVICE 

3  MARKOV  GENERAL 

4  VACANCY 

1 

ENTER  N  (INITIAL  STOCK  VECTOR) 

1  29  74  28  11 

ENTER  THE  PROMOTION  RATE  VECTOR. 

THIS  VECTOR  SHOULD  INCLUDE  THE  PROMOTION  RATES  FOR  THE  FIRST 
3  CLASSES. 

.  102  .  046  .033 

ENTER  THE  WASTAGE  RATE  VECTOR. 

THE  4TH  VALUE  SHOULD  INCLUDE  THE  TOTAL  RATE  FOR  THAT 
STATE  DUE  TO  EITHER  PROMOTION  OR  WASTAGE. 

17  • 1 2 4  • 1  .098 

ENTER  THE  NUMEEE  OF  THE  RECRUIT  TYPE 

1  FIXED  RECRUIT  VECTOR 

2  ADDITIVE  (RECRUIT  SIZE) 

3  MULTIPLICATIVE  (RECRUIT  SIZE) 

4  ADDITIVE  (SYSTEM  SIZE) 

5  MULTIPLICATIVE  (SYSTEM  SIZE) 

4 

ENTER  R  PROP  (RECRUITMENT  PROPORTION  VECTOR) 

10  0  0 

ENTER  ADDITIVE  INCREASE 
0 

ENTER  THE  PERCENT  CODE 
0  NO  GRADE  PERCENTAGES 

1  GRADE  SIZE  AS  PERCENT  OF  TOTAL  GRADE  SIZE 

2  GRADE  SIZE  AS  PERCENT  OF  ORIGINAL  GRADS  SIZE 
7  QUIT  PROGRAM 

1 

WOULD  YOU  LIKE  TO  SEE  THE  ENTERED  DATA? 


0 

0 

0.033 

0.902 


P  MATRIX 

0.728  0.102  0  0 

0  0.83  0.  046  0 

0  0  0.  867  0.033 

0  0  0  0.902 

N  VECTOR 
129  74  28  1  1 
OPTION  =4 
INC  =0 

RECRUITMENT  PROPORTION  VECTOR 
10  0  0 

WOULD  YOU  LIKE  TO  CHANGE  aNY  CF  THE  DATA? 
0  NO 
1  YES 

7  QUIT  PROGRAM 


ENTER  THE  NUMEER  OF  THE  YEAR  YOU  WISH  TO  SEE 
1  0 

DO  YOU  WISH  TO  SEE  THE  INTERVENING  YEARS? 

0  NO 
1  YES 


T 

N 

PERCENT 

0  1 

129 

531 

2 

74 

31) 

3 

28 

12) 

4 

11 

5 

TOTAL 

242 

100) 

10  1 

129 

-  53l 

2 

77 

32) 

3 

27 

11) 

4 

10 

4) 

TOTAL 

242 

100 

DO  YOU  WISH  TC  SEE  ANY  OTHER  YEARS? 
0  NC 
1  Y  ES 

7  QUIT  PROGRAM 


ENTER  THE  NUMBER  OF  THE  YEAR  YOU  WISH  TO  SEE 
50 

DC  YOU  WISH  TO  SEE  THE  INTERVENING  YEARS? 

0  NO 
1  YES 


1 

2 

3 

4 

TOTAL 


DO  YOU  WISH  TO  SEE  ANY  OTHER  YEARS? 
0  NO 
1  YES 

7  QUIT  PROGRAM 


ENTER  THE  NUMBER  OF  THE  YEAR  YOU  WISH  TO  SEE 

100 

DO  YOU  WISH  TO  SEE  THE  INTERVENING  YEARS? 

0  NO 
1  YES 


1 

129 

53 

2 

77 

32 

3 

27 

11 

4 

9  r 

4 

TOTAL 

242  ' 

100 

DO  YOU  WISH  TO  SEE  ANY  OTHER  YEARS? 
0  NO 
1  Y  rS 

7  QUIT  PROGRAM 


DO  YOU  WISH  TO  SEE  THE  STEADY  STATE  VECTOR? 
0  NO 
1  Y  ES 


PERCENTAGES  ARE  IN  STEADY  STATE 

999  1  129  (  53) 

2  .77  I  32 

3  27  (  11 

4  9  (  4 

TOTAL  242  (  100) 


ARE  YOU  THROUGH? 
0  NO 
1  YES 


2.  C^se  J 


Bxampi e  (b) 


START 

DO  YOU  WISH  TO  ENTER  DATA? 

0  NO 

1  YES 

• 

1 

ENTER  THE  NUMEER  OF  THE  MODEL  TYPE 

1  MARKOV  HIERARCHICAL 

2  MARKOV  LENGTH  OF  SERVICE 

3  MARKOV  GENERAL 

4  VACANCY 


ENTER  N  (INITIAL  STOCK  VECTOR ) 

300  200  100 

ENTER  P  (TRANSITION  MATRIX)  El  ROWS 
ENTER  1TH  ROW 

.6  .15  .05 
ENTER  2TH  ROW 

.75  .2  0 
ENTER  3 1H  ROW 

.05  0  .9 

ENTER  THE  NOMEER  OF  THE  RECRUIT  TYPE 

1  FIXED  RECRUIT  VECTOR 

2  ADDITIVE  (RECRUIT  SIZE) 

3  MULTIPLICATIVE  (RECRUIT  SIZE) 

4  ADDITIVE  (SYSTEM  SIZE) 

5  MULTIPLICATIVE  (SYSTEM  SIZE) 

4 

ENTER  REROP  (RECRUITMENT  PROPORTION  VECTOR) 

.7  .07  .23 

ENTER  ADDITIVE  INCREASE 
« 

« 

0 

ENTER  THE  PERCENT  CODE 
0  NO  GRADE  PERCENTAGES 

1  GRADE  SIZE  AS  PERCENT  OF  TOTAL  GRADE  SIZE 

2  GRADE  SIZE  AS  PERCENT  OF  ORIGINAL  GRADE  SIZE 
7  QUIT  PROGRAM 

1 

WOULD  YOU  LIKE  TO  SEE  THE  ENTERED  DATA? 

0  NO 

1  YES 

1 

P  MATRIX 
0.  6  0.15  0.  05 

0.  75  0.2  0 

0.  05  0  0.9 

N  VECTOR 
300  200  100 
OPTION  =4 

INC  =0 

RECRUITMENT  PROPORTION  VECTOR 
0.7  0.07  0.  23 


WOULD  YOU  LIKE  TO  CHANGE  ANY  OF  THE  DATA? 

0  NC 
1  Y  ES 

7  QUIT  PROGRAM 
0 

ENTER  THE  NUMEER  OF  THE  YEAR  YOU  WISH  TO  SEE 
1  0 


69 


DO  YOU  WISH  1C  SEE  THE  INTEBV ENING  YEAES? 
0  NO 
1  YES 


DO  YOU  WISH  TO  SEE  ANY  OTHER  YEARS? 

0  NC 
1  Y  ES 

7  QUIT  PROGRAM 

1 

ENTER  THE  NUMBER  OF  THE  YEAR  YOU  WISH  TO  SEE 
50 

DO  YOU  WISH  TO  SEE  THE  INTERVENING  YEARS? 

0  NO 
1  YES 


50  1 

258  ( 

43) 

2 

54  { 

9 

3 

287  ( 

48) 

TOTAL 

600  ( 

100) 

69 

DO  YOU  WISH  TC  SEE  ANY  OTHER  YEARS? 

0  NO 
1  YES 

7  QUIT  PROGRAM 

1 

ENTER  THE  NUMEER  OF  THE  YEAR  YOU  WISH  TO  SEE 
1  00 

DO  YOU  WISH  TC  SEE  THE  INTERVENING  YEARS? 

0  NO 
1  YES 

0 


100  1 

258 

43 

2 

54  r 

9 

3 

287  i 

48 

TOTAL 

600  ! 

100 

70 


DO  YOU  WISH  TO  SEE  AMI  OTHER  YEARS? 
0  MC 
1  YES 

7  QUIT  PROGRAM 


DO  YOU  WISH  TO  SEE  THE  STEADY  STATE  VECTOR? 
0  NC 
1  YES 

* 

1 


PERCENTAGES  ABE  IN 

STEADY 

STATE 

999  1 

258 

43 

2 

54 

9 

3 

287  r 

48 

TOTAL 

600  i 

100 

ARE  YOU  THROUGH? 
0  MC 
1  YES 

1 


3.  C^sg  2  Example  (a) 


START 

DO  YOU  WISH  TC  ENTER  DATA? 

0  NC 

1  YES 

• 

1 

ENTER  THE  NUMBER  OF  THE  MODEL  TYPE 

1  MARKOV  HIERARCHICAL 

2  MARKOV  LENGTH  OF  SERVICE 

3  MARKOV  GENERAL 

4  VACANCY 

1 

ENTER  N  (INITIAL  STOCK  VECTOR) 

129  74  28  11 


ENTER  THE  PROMOTION  RATE  VECTOR. 

THIS  VECTOR  SHOULD  INCLUDE  THE  PROMOTION  RATES  FOR  THE  FIRST 
3  CLASSES. 

.  102  .046  .033 

ENTER  THE  WASTAGE  RATE  VECTOR. 

THE  4TH  VALUE  SHOULD  INCLUDE  THE  TOTAL  RATE  FOR  THAT 
STATE  DUE  TO  EITHER  PPOMOTION  OR  WASTAGE. 

.17  .124  .1  .098 

ENTER  THE  NUMBER  OF  THE  RECRUIT  TYPE 

1  FIXED  RECRUIT  VECTOR 

2  A  EDIT  I V E  (RECRUIT  SIZE) 

3  MULTIPLICATIVE  (RECRUIT  SIZE) 

4  ADDITIVE  (SYSTEM  SIZE) 

5  MULTIPLICATIVE  (SYSTEM  SIZE) 


ENTER  R PROP  (RECRUITMENT  PROPORTION  VECTOR) 

10  0  0 

ENTER  ADDITIVE  INCREASE 

m 

5 

ENTER  THE  PERCENT  CODE 
0  NO  GRADE  PERCENTAGES 

1  GRADE  SIZE  AS  PERCENT  OF  TOTAL  GRADE  SIZE 

2  GRADE  SIZE  AS  PERCENT  OF  ORIGINAL  GRADE  SIZE 

7  QUIT  PROGRAM 

1 

WOULD  YOU  LIKE  TO  SEE  THE  ENTERED  DATA? 

0  NO 

1  YES 

• 

1 


0.  728  0.102  0  0 

0  0.83  0.  046  0 

0  0  0.  867  0.033 

0  0  0  0.902 

N  VECTOR 
129  74  28  1  1 

OPTION  =4 

INC  =5 


RECRUITMENT  PROPORTION  VECTOR 
10  0  0 

WOULD  YOU  LIKE  TC  CHANGE  ANY  OF  THE  DATA? 

0  NO 
1  YES 

7  QUIT  PROGRAM 

» 

0 

ENTER  THE  NUMBER  OF  THE  YEAR  YOU  KISH  TO  SEE 
1  0 


DO  YOU  WISH  TO  SEE  THE  INTERVENING  YEA  PS? 
0  NO 
1  YES 

0 


T 


N 


PERCENT 


R 


DO  TOO  WISH  TO  SEE  ANY  OTHEB  YEARS? 

0  NC 
1  YES 

7  QOIT  PROGRAM 

1 

ENTER  THE  NUMEER  OF  THE  YEAR  TOO  WISH  TO  SEE 
50 

DO  YOO  WISH  TO  SEE  THE  INTERVENING  YEARS? 

0  NO 
1  YES 


1 

2 

3 

4 

TOTAL 


273 

56 

155 

31 

49 

10 

15 

3 

492  i 

203 

DO  YOO  WISH  TC  SEE  ANY  OTHER  YEARS? 
0  NC 
1  YES 

7  QUIT  PROGRAM 


ENTER  THE  NUMEER  OF  THE  YEAR  YOO  WISH  TO  SEE 
1  00 

°0  Y02  »ISH  TO  SEE  THE  INTERVENING  YEARS? 

0  NO 
1  YES 


1 

2 

3 

4 

TOTAL 


DO  YOO  WISH  TC  SEE  ANY  OTHER  YEARS? 
0  NO 
1  Y  ES 

7  QUIT  PROGRAM 


ENTER  THE  NUMEER  OF  THE  YEAR  YOO  WISH  TO  SEE 
9  00 

DO  YOU  WISH  TC  SEE  THE  INTERVENING  YEARS? 


2537 
1513 
519 
173 
2 


DO  YOU  aiSH  TC  SEE  AMY  OTHER  YEARS? 
0  NC 
1  YES 

7  QUIT  PROGRAM 


ENTER  THE  NUMEER  OF  THE  YEAR  YOU  WISH  TO  SEE 
9  50 

DO  YOU  WISH  TO  SEE  THE  INTERVENING  YEARS? 

0  NO 
1  YES 


1 

2670 

53; 

2 

1593  r 

32 

3 

547 

11 

4 

182  1 2 3 4 

4' 

TOTAL 

4992 

2063] 

WISH  TO 

SEE  ANY  OTHER 

YEARS 

0  NC 
1  YES 

7  QUIT  PROGRAM 


DO  YOU  WISH  TC  SEE  THE  STEADY  STATE  VECTOR? 
0  NC 
1  Y  ES 


PERCENTAGES  ARE  IN  STEADY  STATE 


1 

2 

3 

4 

TOTAL 


2670 

1593 

547 

182 

4992 


ARE  YOU  THROUGH? 
0  NO 
1  YES 


4.  C^se  2 
START 


Example  ( fc) 


DO  YOU  KISH  TC  ENTER  DATA? 
0  NC 


ENTER  THE  NUMEER  OF  THE  MODEL  TYPE 

1  MARKOV  HIERARCHICAL 

2  MARKOV  LENGTH  OF  SERVICE 

3  MARKOV  GENERAL 

4  VACANCY 

3 

ENTER  N  (INITIAL  STOCK  VECTOR) 

300  200  100 

ENTER  P  (TRANSITION  MATRIX)  BY  ROWS 
ENTER  1  IH  RON 

.6  .15  .0  5 
ENTER  2TH  RON 

.75  .2  0 
ENTER  3TH  ROM 

.  05  0  .9 

ENTER  THE  NUMEER  OF  THE  RECRUIT  TYPE 

1  FIXED  RECRUIT  VECTOR 

2  ADDITIVE  {RECRUIT  SIZE) 

3  MULTIPLICATIVE  (RECRUIT  SIZE) 

4  ADDITIVE  (SYSTEM  SIZE) 

5  MULTIPLICATIVE  (SYSTEM  SIZE) 

4 

ENTER  R PROP  (RECRUITMENT  PROPORTION  VECTOR) 
.7  .07  .23 

ENTER  ADDITIVE  INCREASE 
1  25 


ENTER  THE  PERCENT  CODE 
0  NO  GRADE  PERCENTAGES 

1  GRADE  SIZE  AS  PERCENT  OF  TOTAL  GRADE  SIZE 

2  GRADE  SIZE  AS  PERCENT  OF  ORIGINAL  GRADE  SIZE 
7  QUIT  PROGRAM 

1 

WOULD  YOU  LIKE  TO  SEE  THE  ENTERED  DATA? 

0  NO 

1  YES 

1 

P  MATRIX 
0.6  0.15  0.  05 

0.  75  0.2  0 

0.  05  0  0.9 

N  VECTOR 
300  200  100 

OPTION  =4 
INC  =25 

RECRUITMENT  PROPORTION  VECTOR 
0.7  0.07  0.  23 


WOULD  YOU  LIKE  10  CHANGE  ANY  OF  THE  DATA? 
0  NO 
1  YES 

7  QUIT  PROGRAM 


ENTER  THE  NUHEER  OF  THE  YEAR  YOU  WISH  TO  SEE 
1  0 

DO  YOU  WISH  TC  SEE  THE  INTERVENING  YEARS? 

0  NO 
1  YES 


PERCENT 


1 

2 

3 

IOTA  L 

r 

2 

3 

TOTAL 


300 

200 

100 

600  : 

50 

33 

17 

100 

4 77  ~T 

~49 

88 

10 

345 

41 

850  i 

142 

DO  YOU  WISH  TC  SEE  ANY  OTHER  YEARS? 
0  NC 
1  YES 

7  QUIT  PROGRAM 


ENTER  THE  MUMEE8  OF  THE  YEAR  YOU  WISH  TO  SEE 
50 


DO  YOU  WISH  1C  SEE  THE  INTERVENING  YEARS? 
0  NO 
1  YES 


1 

2  • 
3 

TOTAL 


827 

173 

849 

1850 


DO  YOU  rilSH  TC  SEE  ANY  OTHER  YEARS? 
0  NO 
1  YES 

7  QUIT  PROGRAM 


ENTFP  THE  NUMBER  OF  THE  YEAR  YOU  WISH  TO  SEE 
1  00 

DO  YOU  WISH  TC  SEE  THE  INTERVENING  YEAPS? 

0  NO 
1  YES 


v\<v*  ' 


ft 


APPENDIX  E 

COMPUTER  PRINTOUTS  OF  EXAMPLES  FOB  OPTION  5 


1  -  Example  (a) 


S1ART 

DO  IOO  WISH  TO  ENTER  DATA? 

0  NC 

1  YES 

1 

ENTER  THE  NUMBER  OF  THE  MODEL  TYPE 

1  MARKOV  HIERARCHICAL 

2  MARKOV  LENGTH  OF  SERVICE 

3  MARKOV  GENERAL 

4  VACANCY 

1 

ENTER  N  (INITIAL  STOCK  VECTOR) 

129  74  28  11 

ENTER  THE  PROMOTION  RATE  VECTOR. 

THIS  VECTOR  SHOULD  INCLUDE  THE  PROMOTION  RATES  FOR  THE  FIRST 
3  CLASSES. 

.102  .046  .033 

ENTER  THE  WASTAGE  FATE  VECTOR. 

THE  4TH  VALUE  SHOULD  INCLUDE  THE  TOTAL  RATE  FOR  THAT 
STATE  DUE  TO  EITHER  PROMOTION  OR  WASTAGE. 

.  17  .  124  .  1  .098 

ENTER  THE  NUMEER  OF  THE  RECRUIT  TYPE 

1  FIXED  RECRUIT  VECTOR 

2  ADDITIVE  (RECRUIT  SIZE) 

3  MULTIPLICATIVE  (RECRUIT  SIZE) 

4  ADDITIVE  (SYSTEM  SIZE) 

5  MULTIPLICATIVE  (SYSTEM  SIZE) 

5 

ENTER  R PROP  (RECRUITMENT  PROPORTION  VECTOR) 

10  0  0 

ENTER  MULTIPLICATIVE  FACTOR 

1  .01 

ENTER  THE  PERCENT  CODE 
0  NO  GRADE  PERCENTAGES 

1  GRADE  SIZE  AS  PERCENT  OF  TOTAL  GRADE  SIZE 

2  GRADE  SIZE  AS  PERCENT  OF  ORIGINAL  GRADE  SIZE 
7  QUIT  PROGRAM 

1 


79 


WOULD  you  LIKE  TO  SEE  THE  ENTERED  DATA? 
0  NO 


0 

0 

0.033 

0.902 


P  MATRIX 

0.  728  0.102  0  0 

0  0.83  0.  046  0 

0  0  0.  867  0.033 

0  0  0  0.902 

N  VECTOR 
129  74  28  1  1 
OPTION  =5 
FACT  =1.0  1 

RECRUITMENT  PROPORTION  VECTOR 
10  0  0 


WOULD  you  LIKE  TC  CHANGE  ANY  OF  THE  DATA? 
0  NC 
1  YES 

7  QUIT  PROGRAM 


ENTER  THE  NUMBER  OF  THE  YEAR  YOU  WISH  TO  SEE 
1  0 


DO  YOU  WISH  TC  SEE  THE  INTERVENING  YEARS? 
0  NO 
1  YES 


T 

N 

PERCENT  P. 

0  1 

129  | 

[  53) 

2 

74  i 

31) 

3 

28  | 

12) 

4 

11  1 

5) 

TOTAL 

242  | 

:  ioo) 

10  1 

148  l 

t  55) 

2 

82  | 

3 

27  { 

:  io( 

4 

10  i 

4) 

TOTAL 

268  I 

111 

DO  YOU  WISH  TC  SEE  ANY  OTHER  YEARS? 
0  NC 
1  Y  ES 

7  QUIT  PROGRAM 


ENTER  THE  NUMEER  OF  THE  YEAR  YOU  WISH  TO  SEE 
5  0 


CO  YOU  WISH  TO  SEE  THE  INTERVENING  YEARS? 
0  NO 
1  YES 


••  .>  .*•  /• 


50  1 

2 

3 

4 

TOTAL 


222 

55 

125  . 

31 

40 

10 

12 

3 

400 

165 

62 


DO  YOU  WISH  TC  SEE  ANY  OTHER  YEARS? 

0  NC 
1  YES 

7  QUIT  PROGRAM 

1 

ENTER  THE  NUMEER  OF  THE  YEAR  YOU  WISH  TO  SEE 
1  00 

DO  YOU  WISH  TC  SEE  THE  INTERVENING  YEARS? 

0  NO 
1  YES 

0 


100  1 
2 

3 

4 

TOTAL 


362 

55 

206 

31 

66 

10 

20 

3 

654 

270 

101 


DO  YOU  WISH  TO  SEE  ANY  OTHER  YEARS? 
0  NC 

1  V  jl  ? 

1  QUIT  PROGRAM 
0 


DO  YOU  WISH  TO  SEE 
0  NC 
1  Y  ES 

1 

PERCENTAGES  ARE  IN 
909  1 

o 

3 

4 

TOTAL 


THE  STEADY  STATE 


STEADY  STATE 


3  6  2  { 

55) 

2  06  ( 

31) 

DO  ( 

10) 

20  ( 

3) 

654  ( 

100) 

VECTOR? 


ID  1 


ARE  YOU  THROUGH? 
0  NO 
1  YES 

1 


2.  Exajrgle  (b) 


START 

DO  YOU  WISH  TC  ENTER  DATA? 
0  NC 
1  YES 

1 


81 


ENTER  THE  NUMEER  OF  THE  MODEL  TYPE 

1  MARKOV  HIERARCHICAL 

2  MARKOV  LENGTH  OF  SERVICE 

3  MARKOV  GENERAL 

4  VACANCY 


ENTER  N  (INITIAL  STOCK  VECTOR) 

3  00  200  100 

ENTER  P  (TRANSITION  MATRIX)  BY  ROFS 
ENTER  1  TH  ROW 

.6.15  .05 

ENTER  2 TH  ROW 

.75  .2  0 
ENTER  3  TH  ROW 

.  05  0  .9 

ENTER  THE  NUMEEr,  OF  THE  RECRUIT  TYPE 

1  FIXED  RECRUIT  VECTOR 

2  ADDITIVE  (RECRUIT  SIZE) 

3  multiplicative  (recruit  size) 

“  ADDITIVE  (SYSTEM  3IZ*) 

5  TULTI  PJ.lCrfllVE  (SYSTEM  SIZE' 

5 

ENTER  REROP  (RECRUITMENT  PROPORTION  VECTOR) 

.7  .07  .2  3 

ENTER  MULTIPLICATIVE  FACTOR 

1  .01 

ENTER  THE  PERCENT  CODE 
0  NO  GRADE  PERCENTAGES 

1  GRADE  SIZE  AS  PERCENT  OF  TOTAL  GRADE  SIZE 

2  GRADE  SIZE  AS  PERCENT  OF  ORIGINAL  GRADE  SIZE 

7  QUIT  PROGRAM 

1 

WOULD  YOU  LIKE  TO  SEE  THE  ENTERED  DATA? 

0  NO 

1  YES 

1 

P  MATRIX 
0.6  0.15  0.05 

O.  75  0  .2  0 

0.  05  0  0.9 

N  VECTOR 
300  200  100 
OPTION  =5 
FACT  =1.01 

RECRUITMENT  PROPORTION  VECTOR 
0.7  0.07  0.  23 

WOULD  YOU  LIKE  TO  CHANGE  ANY  OF  THE. DATA? 

0  NC 
1  Y  ES 

7  QUIT  PROGRAM 
0 


82 


ENTER  THE  NUMBER  OF  THE  YEAR  YOU  KISH  TO  SEE 


1  0 

DO  YOU  WISH  TC  SEE  THE  INTERVENING  YEARS? 
0  NO 
1  YES 


0 


T 

N 

PERCENT 

R 

. "o' 

1 

300 

50] 

i 

2 

200 

33 

1 

3 

100 

17 

| 

TOTAL 

600 

100] 

1 

10 

1 

3  Id 

r  48] 

I 

2 

67 

10 

1 

3 

277 

42 

| 

TOTAL 

66  1 

1 10 

1 

DO  YOU  WISH  TO  SEE  ANY  OTHER  YEARS? 

0  NC 
1  YES 

7  QUIT  PROGRAM 
1 

ENTER  THE  NUMEER  OF  THE  YEAR  YOU  WISH  TO 
50 

DO  YOU  WISH  TO  SEE  THE  INTERVENING  YEARS? 
0  NO 
1  YES 

0 


50 


1 

2 

3 

TOTAL 


44 

9 

46 

164 


SEE 


123 


DO  YOU  WISH  TC  SEE  ANY  OTHER  YEARS? 

0  NC 
1  YES 

7  QUIT  PROGRAM 
1 

ENTER  THE  N'JMEER  OF  THE  YEAR  YOU  WISH  TO  SEE 
1  00 

DO  YOU  WISH  TO  SEE  THE  INTERVENING  YEARS? 

0  NO 
1  YES 


0 


100  1 

717  ( 

44] 

1 

2 

150  { 

9' 

)  . 

3 

754 

47] 

1 

TOTAL 

1621  ( 

270 

i 

20  3 

8  3 


DO  YOU  WISH  TO  SEE  ANY  OTHER  YEARS? 
0  NO 
1  YES 

7  QUIT  PROGRAM 

« 

0 


DO  YOU  WISH  TC  SEE 
0  NC 
1  YES 

• 

1 

PERCENTAGES  ARE  IN 
999  1 

2 
3 

TOTAL 


THE  STEADY  STATE  VECTOR? 


STEADY  STATE 
717 
150 
754 
1621 


44) 

9) 

47) 

100) 


20 


ARE  YOU  THROUGH? 
0  NC 
1  YES 

1 


84 


REPRODUCED  AT  GOVERUMfU I 


APPEHDIX  P 

1PL  COMPOTES  PHOGBAM  FOR  STEADY  STATE  DISTRIBOTIOM  FOMCTIOH 


This  is  the  APL  function  that  is  used  to  find  the  steady 
state  distribution  vector  for  all  options. 


V^STDST ;  Q ;  IDENT  ,RV-,W ;  NSS;  NHATX  ;  NR  AT  2  ;  FACT  I ;  IQ ;  ZERO 

RV*(X,K)pRPROP 
IDENT*- (K,K)p(X.KpO) 

W*-(X,K)pX-  U>\2)  P) 

Q*-P+(  +  .xRV) 

FAT.lm  UlDENT-  HF1X- ADDREC  ■  K^EEC,ADDSYS ,  MOLTS!  S 
NSS*-N*(+/N) 

N^D  ROOND~N 
+END 

ADDREC :  NHATX*-RPROP+ . x®  ( IDENT- P) 

NSS*-NHATXi (+/NRATX) 

■+END 

KULTREC : FACTI+FACT* IDENT 
NHAT2+R+ . xg (FACTI-P) 

NSS*-NHAT2*{ +/NHAT2 ) 

+END 

ADDSYS : IQ*-IDENT -  Q 


IQ\-,X\+KpX 
ZERO*- 1  .  (K-  1  ^  oO 
NSS*-ZER0  + . x®JQ 


MULT  SYS  :  FA  CTI  *-FACT*I  DENT 
NSS*- (FACT-  1  )  xRPROP+  .  x®  ( FACT  I-  Q  ) 
END  :  TC0UNT*-9  9  9 
PERCENT*- L  (0  .  5  +  lOOx/VSS) 

TOTPERCENT*- 10  0  - 

OUTPUT 

->0 
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